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I   INTRODUCTION 

1.1       The  purpose  of  the  work  reported  in  the  following  is  to  examine  the 

import  of  various  assumptions  on  the  nature  of  two-dimensional  stationary 

hydrodynamic  flow.  We  assume  throughout  that  the  fluid  under  consideration 

is  a  continuum,  that  its  behavior  is  governed  by  the  equations  of  conservation 

of  mass,  momentum  and  energy,  and  that  the  second  law  of  thermodynamics  is 

applicable.  We  further  assume  that  the  flow  takes  place  in  a  plane Tf }      that 

is,  the  flow  variables  are  independent  of  distance  perpendicular  to  TT  and  all 

quantities  are  functions  of  two  coordinates „  We  will  consider  a  section  of 

three  space  parallel  to  ''  which  has  unit  thickness  perpendicular  to  ''  .   In 

the  material  given  below  a  unit  cross  section  of  the  fluid  will  refer  to  an 

element  of  volume  in  this  section,, 

The  conservation  equations  will  be  stated  in  a  Cartesian  coordinate 

system.  Let  V     be  a  closed  Jordan  curve  fixed  with  respect  to  the  set  of 

coordinates  x  and  x  .   Specifically  the  conservation  equations  state  that  the 

total  rate  of  mass,  momentum  or  energy  entering  the  region  bounded  by  P  is 

zero,  the  force  acting  on  the  boundary  or  the  work  being  done  on  the  boundary 

respectively.  The  rate  of  transfer  of  mass  across  a  unit  straight  line  whose 

unit  normal  has  direction  cosines  n.  is  pu.a.  so  that  the  conservation  of  mass 

1    t-  1  x 


is 


f 


^u.n.ds  =  0  (1.1.1) 

V 

where  p  is  the  mass  per  unit  volume  of  an  x  -  x  unit  cross  section  of  the  fluid, 
u.  are  the  components  of  the  velocity  referred  to  the  Cartesian  coordinates, 
and  the  summation  convention  is  used  for  Roman  suffices  running  through  1 
and  2. 


The  ith  component  of  momentum  flow  across  a  unit  line  is  pu .u.n.  and  the  ith 


component  of  the  force  acting  on  a  unit  line  is  T   .n.  where 


so  that  the  conservation  of  momentum  is 


n 


where     w  is  the  viscosity  of  the  gas, 


p  the  pressure 

and  Y   .   means     ^ — 

,J  dx. 


/<fVj  +  1Ws  =  °  (1-1-2) 


The 


2 


rate  of  flow  of  energy  across  the  unit  line  is  pn.u.v  /2  due  to  the  kinetic 


S 


energy  and  pu.n./(^-  l)  due  to  the  internal  energy,  so  that  the  conservation 
of  energy  can  be  written  as 

[(r^  +  _^)  r.+  T.A   n.u.ds  =  0  (1.1.3) 

LV(  2    y~  l  '   °  ij  ij   j  i 

r 

where     "6      is  the  ratio  of  specific  heats  of  the  gas, 
v  is  the  velocity  of  the  gas 

and  the  heat  conductivity  is  assumed  zero. 

Ji 

The  second  law  of  thermodynamics  implies  that  the  entropy  p/9    will  never 

decrease  in  an  area  moving  with  the  flow.  Throughout  we  will  assume  that  the 
flow  variables  o   ,  u.  and  p  are  bounded. 

1.2       First  we  will  make  the  additional  assumptions  that  the  viscosity  w  is 
zero  and  that  the  velocity  and  pressure  are  single-valued  except  on  a  finite 
number  of  piecewise  analytic  curves.   If  AB  is  an  analytic  segment  of  one  of 
these  curves  with  direction  cosines  \.}    and  if  PQ  and  SR  are  representations 
of  AB  displaced  a  distance  d  along  the  normal  to  AB  at  each  point,  to  the 


Figure  1,1 

left  for  PQ  and  the  right  for  SR  such  that  PQRS  encloses  no  other  dis- 
continuities,, then  we  can  write  the  conservation  equations  in  the  form 


Xds=       /      Xds+      j      Xds+  Xds+       J      X  ds     =     0 

SR  QP  RQ  PS 


where  X  can  be  the  mass  flow,  momentum  flow  plus  force,  or  the  energy  flow  plus 

work.   Since  this  is  bounded  we  have  as  d  tends  to  zero   I   X  ds  =  J  X  ds, 

"SR        PQ 
or  writing  X  for  the  limit  on  the  left  of  AB  and  X   for  the  limit  on  the 

right  we  get 


/   (X+  -  X")  ds  =  0 
J  A 


Since  this  is  true  for  any  points  A  and  B  such  that  AB  is  analytic,  we 
have  X  =  X   almost  everywhere.   These  are  the  Rankine-Hugoniot  equations,  or 
R-H  equations,  and  imply  that  either  the  discontinuity  is  a  stream  line  of  the 
flow  and   [p]  =  0  almost  everywhere  where    [p]  =  p  -  p  ,  or  that  the  dis- 
continuity is  a  shock  described  by  the  equations 


ut  -  u;  d„2a) 

_2 

r 


^    =    y  =  i  +  yfi  (<-2  -  i)  (1.2.2) 


A  =  f  =      LT^lgl  (1.2.3) 

f"  <■  2  +  (/-l)cr2 


k 

u 

<r    ^ 

almost  everywhere,  where  c   =  —  ,   cr=  —  >   ^  an(^  ^  are  the  components 

r  c 

of  the  velocity  tangential  and  normal  to  the  shock  respectively,  and  the  super- 
fix  has  been  dropped.  The  second  law  of  thermodynamics  implies  that  y  >  1, 
or  equivalently  that  o~  >   1. 

In  regions  where  the  flow  variables  are  different iable  we  will  transform 
the  independent  variables  x  and  x  to  s  and  T  where  T  is  defined  as  0  along  an 
arbitrary  curve  which  is  nowhere  tangential  to  a  stream  line  and  meets  each 
stream  line  once  and  only  once.  The  T  coordinate  of  any  point  is  the  arc  length 

from  T  =  0  along  the  unique  stream  line  to  that  point  so  that  the  transforma- 

^x.   u. 

tion  is  given  by  rrr-  =  —  and  the  value  of  x  (s,  0)  and  x  (s,0)  as  functions 
o  T     V  x  c. 

of  s.  Thus  the  s  coordinate  is  defined  for  a  stream  line  and  hence  for  each 
point  on  it  by  prescribing  s  along  the  curve  T  =  0.  We  define 

Px. 


\ 



l 

1 

Qs 

x2 

= 

q/  -  hh 

u 

= 

JV* 

sin  (j. 

= 

i  i 
\v 

cos  u 

= 

u.  6  .. 
J         *.v 

h 

= 

Y   logp 

r   -  —  to  be  the  angle  between  a  stream  line  and  the  x  axis 

m  =(£"„)  '   =  —  to  be  the  local  Mach  number 
V  3       c 

and      c  =    /—  is  the  local  velocity  of  sound. 


Making  these  substitutions  into  the  conservation  equations  (l.l.l)  -  (1.1.3) 
Taub     has  shown  them  to  become 


where 


7i 


6>s 


=     a.  . 


3h 


ij      Qs 


i  =  1,   2 


sin  u 

2          ~| 
-m  cos  [i 

a.  . 

=     \ 

(—  -  i)cos  n 

sin  u 

m 

(1.2.5) 


(1.2.6) 


dm 
c?T 


1  /--.     r- 1  2v      )i 


m 


(1.2.7) 


Y  -  1 


9T 


(1.2.8) 


plus  the  two  geometrical  relationships 


and 


9\  % 

"5*  =    cos  »    "5T 

_£u_  g?2  sin  u      <?J2 

e>T  9T     "  \  9s 


(1.2.9) 
(1.2.10) 


P  p       p 

When     det  I  |a.  .  |  |   =  X.   (l  -  m  cos  \x)     is  non-zero   (as   it   certainly  is 

when  m  <  l)      | |a. . | |    can  be   inverted  to  give 


^  . 


2T 


fll 

vj       ps 


A      . 


V=  1,    2,    ...,    6,  (1.2.11) 


'J  \(l-m2cos2u) 


sm  u 

1 


1  -  —  cos   n 
m 


2 
m  cos  p. 


sm  a 


-[2  +   (r-l)m2]sin  (i  m2[2   fr   (^-l)m2]   cos   u 


(l -   )   cos   U 

m 

0 


-  —  sm  Li 
m 


2      .  2 

-m     sm  |i  cos     u 


cos    Li. 

-mc   cos  u 
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There  exists  a  class  of  time  dependent  flows  in  two  dimensions  where 

the  time  t  enters  only  in  the  variables  a  =  x  /t  and  a  =  x  /t.  This  is 

called  pseudo- stationary  flow.   If  a  pseudo-stationary  velocity  vector  U.  is 

defined  as  u.  -  a.,   then  if  U.  and  a.  are  used  as  variables  instead  of  u.  and 
I    i  -  i      i  1 

[11] 
x.  the  equations  of  flow  have  been  shown  by  Taub    to  be  of  the  above  form 

with  additional  terms  not  involving  derivatives,  so  that  they  can  be  written  as 

%     -     *„  ^+  =V 

2T     yj  ^s    v 

where  Hy   is  a  function  of  the  flow  variables.  At  the  end  of  Chapter  II  the 
analogy  between  stationary  and  pseudo-stationary  flow  will  be  discussed. 
1.3      These  stationary  and  pseudo- stationary  equations  are  elliptic  or 
hyperbolic  as  the  Mach  number  m  is  smaller  or  larger  than  1.  For  a  Mach 
number  of  exactly  1  they  are  parabolic. 

The  boundaries  of  the  flow  are  combinations  of  shocks  and  walls  which 
do  not  in  general  completely  surround  the  region  being  considered.  However  the 
problem  can  still  be  well  defined  if  the  subsonic  region  can  be  bounded  by  a 
combination  of  shocks,  walls  and  the  sonic  line.  The  supersonic  region  does 
not  need  a  complete  boundary  since  the  equations  are  hyperbolic  there  so  that 
an  initial  value  problem  can  be  solved. 

Typical  configurations  that  have  been  observed  in  shock  tubes  and 
wind  tunnels  are  the  regular  and  Mach  reflection  patterns  and  the  shock  attached 
to  or  detached  from  a  wedge.  In  reflection  an  incident  shock  I  is  reflected 
at  a  wall  W  in  one  of  two  ways.  For  some  incident  shock  angles  a  second  shock 
R  can  be  found  that  will  turn  the  flow  between  I  and  R  parallel  to  the  wall. 
(Figure  1.2)  For  other  incident  angles  the  pattern  of  Figure  1.3  is  observed. 
Here  the  junction  T  of  the  incident  and  reflected  shock  has  moved  away  from 


W      ttt})}  U  /Y TTTTTTTTTTTTTl  7 


Figure  1.2 


w    j  //////  y  /////////  / 

Figure  1.3 
the  wall  and  is  connected  to  the  wall  by  the  Mach  shock  M.  This  is  the  Mach 
reflection  pattern. 

The  problem  of  a  wedge  in  a  uniform  flow  is  similar  in  may  ways.   For 
some  oncoming  flows  and  wedge  angles  attached  shocks  AB  are  possible.   (See 
Figure  1.^)  This  corresponds  to  the  existence  of  the  regular  reflection  of 


Figure  l.k 


8 
Figure  1.2.  The  oncoming  flow  in  the  wedge  problem  corresponds  to  the  flow 
behind  the  incident  shock  in  the  regular  reflection.  The  function  of  the 
attached  shock  or  the  reflected  shock  is  to  turn  the  flow  through  an  angle 
so  that  it  is  parallel  to  the  wall  or  the  wedge.   If  we  think  of  the  wall  as 
being  the  side  AC  of  the  wedge,  then  the  problems  are  identical. 

When  an  attached  shock  is  no  longer  possible  a  detached  shock  can  be 
observed  in  front  of  the  wedge  that  must  result  in  the  velocity  being  zero  at 
the  apex  of  the  wedge.   This  corresponds  loosely  to  a  three  shock  configuration, 
the  Mach  reflection  phenomenum,  where  the  triple  point  T  detaches  from  the 
wall,   (See  Figure  1.3)«  In  the  three  shock  configuration  there  is  always  a 
trivial  configuration  where  the  reflected  shock  R  is  a  Mach  wave  and  the  Mach 
shock  M  is  a  continuation  of  the  incident  shock  I.   If  the  incident  shock  is 
sufficiently  weak  this  is  the  only  theoretical  solution,  but  experimentally 
other  solutions  have  been  observed.  We  shall  return  to  this  point  in  subsequent 
discussions. 

In  Chapter  II  we  will  use  the  equations  of  Section  1.2  to  study  the 

behavior  of  non- viscous  flow  around  a  wedge  in  a  supersonic  oncoming  flow. 

[2] 
This  is  equivalent  to  the  problem  of  regular  reflection.   Buseman    has 

surveyed  the  results  that  have  been  obtained  in  the  wedge  problem.  When  the 

wedge  angle  is  small  there  are  two  shocks  possible,  one  weak  and  one  strong. 

He  points  out  that  both  have  been  obtained  experimentally.  As  the  wedge  angle 

&  increases  the  angles  ccn  of  these  shocks  approach  each  other  until  they  are 

equal  at  the  so-called  extreme  angle.   For  larger  %    no  attached  shock  is 

possible. 

If  %    is  small  enough  the  flow  behind  the  weak  shock  is  supersonic 

and  the  shock  is  straight  from  the  apex  A  of  the  wedge  (see  Figure  1„4)  to  the 

point  B  where  the  characteristic  leaving  the  shoulder  C  of  the  wedge  meets 


9 
the  shock.  From  there  it  consists  of  a  number  of  analytic  segments  at  least 
until  the  flow  becomes  sonic.   In  Section  2.1  this  will  be  shown  in  some  detail. 
As  &  increases  the  angle  of  the  weak  shock  an  increases  until  the  flow  behind 
the  shock  is  sonic  at  a  =  a   the  sonic  angle.   If  %    is  increased  still 

\J  s 

[12] 
further  the  shock  is  curved  all  of  the  way  down  to  the  point  A.  Thomas     and 

Taub    point  out  that  only  for  a  countable  number  of  angles  "Xs    is  it  possible 

to  have  an  analytic  shock  that  is  curved  and  yet  have  a  straight  stream  line  AC. 

Guderly       has  studied  the  singular  flow  in  ABC  when  the  oncoming  flow  is 

nearly  sonic  and  has  obtained  a  solution  such  that  the  curvature  of  the  shock 

M 

at  A  is  zero  for  a_  between  the  sonic  angle  and  the  Crocco    angle.  At 
the  Crocco  angle  the  shock  may  have  a  finite  curvature  and  the  solution  be 
regular.  For  a_  greater  than  this  but  less  than  the  extreme  angle  which  is 
the  a       corresponding  to  the  largest  possible  deflection  ~)C  of  the  flow,  the 
shock  may  have  an  infinite  curvature  at  A. 

Later  in  Chapter  II  we  discard  the  assumption  that  the  shock  is  piece- 
wise  analytic  when  it  bounds  a  subsonic  region,  and  introduce  a  non-analytic 
transformation  of  the  s-T  plane.   In  this  new  plane  we  assume  that  we  have  a 
singularity  that  has  the  nature  of  a  branch  point.  The  effect  of  this  assumption 
is  that  in  the  physical  plane  the  flow  variables  are  assumed  single-valued  behind 
the  shock  although  their  derivatives  have  values  at  the  singularity  that  depend 
on  the  direction  in  which  it  is  approached.  The  angle  of  the  shock  thus  has  a 
limit  as  one  approaches  the  apex  A  of  the  wedge  and  this  angle  must  be  such 
that  the  flow  emerging  from  the  shock  is  parallel  to  the  wedge  wall  at  A.  This 
does  not  allow  the  type  of  singularity  where  the  pressure  or  direction  of  flow 
is  multiple -valued  behind  the  shock  to  exist .   Consequently  the  angle  of  the 
shock  at  the  apex  of  the  wedge  is  determined  by  the  R-H  equations  and  the  fact 
that  the  flow  leaving  the  shock  must  be  parallel  to  the  wall. 
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Singularities  of  the  type  treated  in  Chapter  II  cannot  be  used  in 
discussing  three  shock  configurations  occurring  in  Mach  reflection  of  weak 
incident  shocks  since  one  would  have  to  postulate  the  existence  of  a  non-trivial 
three  shock  configuration  in  order  to  place  such  a  singularity  at  the  triple 
point  of  the  configuration  of  shocks  in  this  case. 

The  local  boundary  conditions,  that  is  the  flow  ahead  of  the  shock 
and  the  shape  of  the  wall  in  the  neighborhood  of  the  apex  of  the  wedge,  determine 
a  class  of  solutions  in  the  new  plane.  This  class  includes  the  regular  solutions 
and  others.   The  solutions  to  some  of  the  boundary  value  problems  may  require 
shocks  having  an  infinite  curvature. 

1.^       One  problem  that  this  method  will  fail  to  provide  a  solution  for  is 
that  of  a  curved  wedge  in  a  uniform  oncoming  flow  when  the  angle  of  the  wedge 
satisfies  a  certain  condition.  There  are  a  countable  number  of  such  angles  and 
unless  additional  restrictions  are  satisfied  by  the  flow  and  wedge  profile  there 
are  no  solutions  of  the  class  discussed.   Rejecting  the  idea  that  nature  will 
discriminate  against  this  set  of  wedge  angles  we  must  conclude  either  that  there 
are  solutions  with  other  forms  of  singularities  or  that  the  mathematical  model 
of  the  fluid  is  inadequate. 

One  type  of  singularity  that  we  might  look  for  is  one  where  the  angle 
of  the  flow  and  other  flow  variables  are  not  single- valued  everywhere  behind 
the  shock.   Such  a  flow,  the  Prandtl-Meyer    flow,  exists  in  the  supersonic 
region,  but  we  are  concerned  with  the  subsonic  region  where  no  such  flow  is 
known  at  present.  We  can  try  and  find  such  singularities  and  try  and  decide 
if  the  mathematical  model  is  adequate  jointly  by  examining  viscous  flow  and 
in  particular  studying  the  solution  as  the  viscosity  tends  to  zero. 

There  is  a  reason  to  believe  that  the  solution  of  a  problem  for  a 
finite  viscosity  will  not  approach  the  solution  of  the  equivalent  non-viscous 
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problem  as  the  viscosity  tends  to  zero.  This  is  that  the  curvature  of  the  shock 
in  some  of  the  solutions  discussed  can  be  infinite.   If  the  flow  is  one-dimen- 
sional, that  is  if  all  derivatives  tangential  to  the  shock  are  zero,  then  the 
R-H  equations  are  nearly  satisfied  across  a  shock  region,  where  by  a  shock 
region  we  mean  the  region  between  two  bounding  lines  T  =  T_  and  T  =  T   which 
encompass  most  of  the  change  in  the  flow.  X  -  X  =  f(T,  -  T_)  where 
f(T..  -  T  )  tends  to  zero  as  T,  -  T_  tends  to  infinity.   If  we  assign  a  value 
to  f  (T.  -  T_)  then  T..  -  T  ,  the  shock  thickness,  is  proportional  to  the  viscosity 
w.  For  small  shock  curvatures  k  terms  will  be  present  that  are  proportional  to 
the  curvature.  The  region  over  which  the  terms  act  is  proportional  to  w  so  the 
effect  is  proportional  to  kw.  Thus  if  the  viscosity  tends  to  zero  the  behavior 
the  solution  will  depend  partly  on  the  behavior  of  kw,  and  if  k  becomes  infinite 

when  the  viscosity  is  zero,  the  terms  in  kw  may  still  play  a  part. 

[10] 
Sternberg    considered  the  effects  of  modifying  the  R-H  equations  in 

the  neighborhood  of  the  intersection  region  by  introducing  viscosity.  Using 

approximate  solutions  obtained  from  an  electric  tank  he  showed  that  a  large 

modification  was  necessary  to  have  an  appreciable  effect  on  the  flow. 

In  Chapter  III  the  equations  of  viscous  flow  will  be  discussed.  To 
simplify  the  problem  we  will  assume  that  in  some  regions  viscosity  can  be 
neglected,  and  that  these  regions  can  be  joined  together  by  transition  regions 
in  which  viscous  forces  play  a  role.  We  will  start  by  examining  the  transition 
region  analogous  to  a  shock,  and  then  examine  ways  in  which  two  or  three  of 
these  regions  can  interact. 

The  treatment  will  differ  from  Sternberg's    treatment  in  that  the 
Mach  and  reflected  shocks  will  be  allowed  to  be  curved  rather  than  being  con- 
sidered as  one- dimensional  shocks.  The  intersection  region  where  all  of  these 
shocks  interact  will  correspond  to  his  non-R-H  region.  He  applied  the  four 
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conservation  equations  around  this  region  to  evaluate  its  size  and  shape  using 
experimental  values  for  the  angles  of  the  three  shocks,  whereas  this  treatment 
will  attempt  to  determine  the  angles  of  the  shocks . 

Sternberg  discussed  the  possible  effects  of  this  region  on  the  flow 
downstream  and  in  his  specific  calculation  showed  that  the  region  is  of  con- 
siderable size  compared  to  the  thicknesses  of  the  shocks.  He  further  argued 
that  the  limit  of  the  solution  as  the  viscosity  tends  to  zero  may  not  be  the 
solution  for  zero  viscosity.  We  too  will  find  that  this  can  be  the  case  for 
regular  and  Mach  reflection. 

Since  the  arguments  presented  are  local  in  character,  no  account  is 
taken  of  distant  boundary  conditions .  The  lengths  involved  are  proportional 
to  viscosity  and  the  shock  curvatures  are  inversely  proportional  to  viscosity. 
By  using  polynomial  approximations  in  s  and  T  for  the  flow  variables  we  will 
study  the  problem  in  detail  and  exhibit  the  equations  governing  two  and  three 
shock  configurations  for  a  linear  polynomial  approximation  to  each  of  the  flow 
variables  p,  Ej  ,  v,  \  and  a.  The  unknowns,  which  are  now  the  coefficients 
of  the  polynomials,  will  be  expressed  as  power  series  in  the  strength  of  the 
incident  shock  so  that  the  equations  for  regular  and  Mach  reflection  can  be 
expressed  as  algebraic  equations  in  the  coefficients  of  the  power  series. 

For  weak  shocks  we  can  take  the  first  terms  of  these  power  series  to 
get  relations  involving  eight  unknowns.  Four  of  these  will  be  variables 
specifying  the  angles  of  the  various  shock  boundaries.  The  solution  for  these 
is  independent  of  the  viscosity  w.  The  other  four  variables  are  two  shock 
curvatures  and  two  distances.   If  we  write  K  =  kw  and  d  =  D/w  where  D  is  a 
typical  length  we  will  find  that  the  solutions  for  K  and  d  are  independent  of 
viscosity  so  that  if  K  /  0,  the  limit  of  the  solutions  as  viscosity  tends  to 
zero  will  be  different  from  the  solution  of  the  non-viscous  shock  intersection 
problem. 
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II  NON- VISCOUS  FLOW 

2.1       In  this  section  the  solution  of  the  stationary  hydrodynamic  equations 
are  discussed  for  a  wedge  in  a  uniform  oncoming  flow  which  is  such  that  the  flow 
behind  the  shock  is  initially  supersonic.   If  the  wedge  has  infinite  straight 
sides  then  a  solution  of  this  type  is  known  to  exist,  namely  the  straight  weak 
shock  that  turns  the  flow  through  the  angle  of  the  wedge.   If  now  a  convex  angle 
appears  in  the  boundary  at  a  point  C,  such  that  the  wall  is  composed  of  two 
straight  portions,  the  shock  will  remain  straight  until  a  characteristic  leaving 
the  point  C  reaches  it  at  B. 

C, 

B 


Figure  2.1 

The  flow  in  the  region  ABC  is  uniform.  Let  BC..D  be  another  characteristic  leaving 
the  point  B.   Then  the  flow  in  the  region  CBO.D  can  only  be  affected  by  the  on- 
coming parallel  flow  from  ABC  and  the  wall  between  C  and  D.   In  the  case  under 
consideration  a  solution  is  known  for  BCD,  namely  a  suitable  Prandtl-Meyer 
region  BC  C  and  a  uniform  region  CO.D.  Hence  all  derivatives  of  the  flow 
variables  taken  along  the  characteristic  BD  can  be  computed.  Now  it  will  be 
shown  that  the  equations  of  motion  and  the  R-H  equations  have  a  solution  that  is 
separately  analytic  in  the  regions  C  BC  and  CpBC  ,  where  BC  is  the  stream  line 
through  the  point  B. 


Ik 


The  coordinates  s  and  T  are  chosen  in  the  region  C  BC.  so  that  T  =  0 

is  the  line  BC  .   In  the  physical  plane  the  characteristics  subtend  an  angle  M 

with  the  stream  line,  where  M  =  sin  (l/m).   Thus  on  BC   u  =  M  -  "#/2  where 

u  is  the  angle  between  the  normal  to  T  equal  constant  (the  curve  BC, )  and  the 

tangent  to  s  equal  constant  (the  stream  line).   Since  derivatives  are  known 

along  BC^  the   - —  are  known  there.   By  (1.2.5)  we  have 
l      <9S 


3s 


=  a 


Hi 

ij  9T 


(2.1.1) 


where 


-cos  M 

1 

sin  M 

a.  . 

=     \ 

2 
-cos  Msin  M 

-   cos  M 



(2.1.2) 


as^follows  from  (1.2.6)  on  substituting  for  u  and  m.  Det  j |a. . | |  =0,  therefore 
>2 


ao  J.  V. 

3s 


and 


9s 


evaluated  along  BC,  will  satisfy  the  relation 


Hi 

9s 


11 


a 


12 


3s 


a. 
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*22 


sin  M  cos  M 


Thus  when  we  have  evaluated  either  -— —  or  -~ —  from  the  Prandtl- Meyer  region 

BCC.  we  can  solve  (2.1.1)  for  -rr=-  •   Since   I  la.  .  I  is  of  rank  1  there  will  be 
1  3T  '  '  ij  ' ' 

a  one  parameter  family  of  solutions.  Therefore  we  can  write 


X 


<Pt     j   n.  j 


(2.1.3) 
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where 


3A 

3s 


\   cos  M 


and 


cos  M 


sin  M 


for  any  constant  (3.. .   Since  the  region  adjacent  to  BC,  is  not  in  general  uni- 
form,  -5 —  £     Oo  Between  the  two  regions  of  analytic  flow  separated  by  the 


<?s 


stream  line  s  =  0  the  flow  variables  F    and  7    must  be  such  that 
(   (T,0)  =   F  (T^O)  for  all  points  on  the  stream  line.   In  particular 


2?n 


evaluated  in  one  region  must  be  the  same  as  that  evaluated  in  the  other  region 
at  B,  for  all  n  =  0,  1,  2,  „ .  „   ■>  For  n  =  0,  this  implies  that  the  initial 


angle  of  the  curved  shock  starting  from  B  is  the  same  as  that  of  the  straight 


shock.  For  region  C  BC.  we  define  T  =  0  as  the  shock  BCL.  Following  Taub 


111] 


we  have  by  differentiation  of  the  R-H  equations 


~    =  ka. 
ds  1 


(2.1.4) 


where  k  is  the  curvature  of  the  shock  and 


2m  sin2u 

*i   ~     r+  1 


a^ 


4  2 
sin  u  /  4 


B*(^T*1?).  eo^-l) 


Inverting   (1.2. 5)  and  substituting  from  (2.1o4)  we  get 

3h 


—     i£A.  . a . 


where 


A,(l  -  m  cos  [i) 


sm  u 


(1- 


cos  u 


m 


(2.1.5) 


2 
m  cos  a 


sin  u 
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Equating  (2.1.3)  and  (2.1.5)  we  get 


kA,  . a .   =  b .  -  B_  c  _, 

ij  J      i    1  i 


(2.1.6) 


These  are  two  equations  in  the  two  unknowns  k  and  f3  .   They  have  a  unique 


solution  if  and  only  if 


det 


A.,  .a. 


A_.a. 
2j  J 


t     0  . 


Since  0  <  u  <  7^2  in  region  BC  C  ,    each  entry  in  A. .  is  positive.  Further  a 

and  a„  are  hoth  positive  in  the  supersonic  region  so  that  A^  .a.  and  A^.a.  are 
2         *  *  &  ±J   J      2j  j 

both  positive.   Since  c  and  c  have  opposite  signs  (2.1.7)  will  always  be 
satisfied. 

Therefore  (2.1.6)  can  always  be  solved  uniquely  to  give 


k  = 


blC2  "  Clb2 


A. .a.c0  -  A_  .a.c. 
lj  j  2   ^2j  j  1 

b1  A_  .a.  -  b0A-.  .a. 
r   =   1  2j  j    2  lj  j 

lj  j  2    2j  j  1 

Now  suppose  that  all  the  partial  derivatives  of  the  flow  variables  of 
order  less  than  n  are  known,  together  with  derivatives  of  k  of  order  less  than 
n-1  such  that  the  equations  connecting  the  derivatives  of  k  and  the  flow  variables 
are  consistent.   By  successive  differentiation  of  (1.2.5)  and  (1.2.6)  we  find 
that 


i1  J.         sn> . 


2s 


n 


1J 


<PT 


n 


(2.1.8) 


where  the  Q  .  are  polynomials  in  the  known  lower 


order  derivatives 


5^i 
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At  B  in  region  C.BCL      the 

1  2       9sn 


axe  known,  therefore 


=  =  P  c.  -  b  . 

9Tn  n  X         ni 


(2.1.9) 


b  .  = 

ni 


9nj 

ir-Sa. 


(-2cos  M)n 

0 
We  get  by  differentiation  of  the  Rankine-Hugoniot  equations  along  the  shock  BC, 


=  a.  =- 

„  n       i    n-1 


<?s 


+  P  . 
l   ~  n-1      ni 


(2.1.10) 


where  P  .  is  a  polynomial  in  the  known  derivatives 
of  lower  order. 
In  the  region  BC C    | |a. . | |  is  not  a  singular  so  (2.1.8)  can  be  inverted  to 


9^ 


=     A' 


% 


'•% 


5  s11 


1J 


(2.1.11) 


,n     ~n-l_ 
=  A. .a.  d      k 

°  °  ~   n-1     ni 

<>^S 


Equating  (2.1.9)  and  (2.1.11)  we  again  have  two  equations  in  two  unknowns  =- 

/-*   n—  -L 

and  3  .  These  have  a  unique  solution  if  and  only  if 


det 


A11 

A.  .a. 


A11 

A.  .a. 


^  0   . 


(2.1.12) 


Following  Taub    we  write 


tanhQ 


,  2   ,  xl/2    . 
=  (m  -  1;    cot  u  , 


hence 


An 


ij  2n/_  2       2  7n/2 

°  \   (1  -  m  cos  uj  ' 


18 


2 
cosh  n9  -   .,  i     sinh  n© 


d-»V^ 


(m  -1) 


m 


sinh  ji9  cosh  in© 


In  the  coordinate  system  lender  consideration 

2        2  77" 

m  cos  \x     <     1     f  0     <     u     <    — —     , 

hence  each  entry  in  A. .  is  positive „  Also 

-*-J 

2  .  2 
_  m  sin  \i 

al  "  y+  1         >     U  ' 
and  a_  may  be  rewritten  as 

v  ,  i   t(l  -  m  cos  |i)  (1 —  +  — —  )  +  (m  -  1)  cos  u  (l  +  — ^-)]  >  ( 

2 
when  m  >  1.   Since  c.  and  c  have  opposite  signs,  (2.1.12)  cannot  vanish  in 

the  supersonic  region  for  n  =  1,  2,  3>  »..   .  Thus  (2.1.9)  and  (2.1.11)  can 

always  be  solved  uniquely  for  — =-  and  P  ,   completing  one  step  of  the 

9s 
induction  process. 

The  region  behind  the  shock  may  be  systematically  covered  by  regions 

in  which  the  flow  is  analytic.  These  regions  start  either  at  the  shock  (i.e. 

regions  like  BCLCL)  or  at  a  characteristic  leaving  a  discontinuity  (i.e.  regions 

like  BC.  C  ) .  They  are  bounded  by  a  stream  line  (which  may  be  the  wall )  and 

terminated  by  a  characteristic  of  the  other  set  leaving  a  discontinuity.   By 

equating  conditions  across  the  stream  line  joining  adjacent  regions  the  flow 

may  be  determined  analytically  in  these  regions.  The  sides  of  the  wedge  can 

consist  of  a  number  of  curved  analytic  segments  instead  of  the  straight  sides 
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treated  earlier,  and  until  another  shock  forms,  or  the  characteristics  of  the  same 

set  intersect  again  (they  intersect  at  C  in  the  Prandtl-Meyer  fan),  the  problems 

can  be  treated  by  the  methods  given. 

2.2       When  the  flow  behind  the  shock  is  subsonic  we  will  have 

det  |  ja.  .  |   =  X,  (l  -  m  cos  u)  >  0 

so  that  |  |a. . J  j  can  be  inverted.  Multiplying  (2.1.8)  by   |  ja. . J  |   =  A. .  we 


get 


*"£   .»       *nJi 


—*•     AT.   =  —  +  P  , 

where  the  P  .  are  polynomials  in  the  known  lower  order  derivatives.   Substituting 


from  (2.1.10()  into  this  we  get 


_±L  =  A«  a   3_Jk  +  An     .  p  (2.2J) 

[12 1        [ill                                         ^  -?2 
Thomas    and  Taub    have  used  this  equation  and  the  fact  that  the  

3Tn 

are  given  at  s  =  0  for  n  >  1  along  a  wedge  wall  or  known  stream  line  to 


show  that  unless  A-„a.  ^  0  for  all  n  >  1,  no  finite may  exist  to 

2j  j  r  _  ^i 

satisfy  (2.2.1)  for  all  n. 

If  the  known  stream  line  is  straight  we  have  that  =  0  for  all 

n 
o>T 

n  >  1.   If  we  assume  that  the  first  n-1  derivatives  of  all  flow  variables  and 

the  first  n-2  derivatives  of  k  are  zero  then  the  P  .  and  the  P  .  of  (2.2.1), are 

nj  nj 

zero,    so  that    (2.2.1)  reduces  to 


r\  -n-: 


An  a       §     y  +     An  a       ^ -       -     0 

A2lai  .  n-1  +  A22a2  n  n-1   "  U  ' 
ds  ds 

n  dn~\. 

Unless  A  a.  =  0  we  must  have  r  =  0  or  00  „   Consequently  we  can  show  by 

J  J  as11"1 

induction  that  unless  A^,. a.  =  0  for  some  n  then  all  derivatives  of  the  curvature 

2j  J 
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[12] 
and  the  curvature  itself  are  zero  or  infinity.  Thomas    has  shown  that  the 


number  of  wedge  angles  for  which  this  happens  is  countable,  and  that  these  wedge 
angles  form  a  dense  set. 

To  study  the  flow  in  the  neighborhood  of  the  apex  A  of  the  wedge  we  will 
make  a  singular  transformation  of  the  s-T  plane.   Let  m,  be  the  Mach  number  of 
the  flow  behind  the  shock  at  A;  that  is  it  is  the  limit  of  the  Mach  number  of 
the  flow  as  one  approaches  A  along  the  shock.  Let  p.,  be  the  value  of  u  at  that 


point . 
Define 


X  =  -f-  -  ■  ^ 


p      p      p     P  p 

R   =  s  (l  -  hl  )  sin  0  +  (s  cos  0  +  T) 

,,    2.1/2  .  u 
s(l  -  m  )  '  sin  p 

tan  (3  =  


s  cos 


0  +  T 


Since  the  shock  cannot  be  a  stream  line  we  know  that  sin  0  /  0  and  we  are 

2 
assuming  that  m   <  1„ 

If  we  introduce  a  further  coordinate  plane  w-z  where 

.      jtft  2Nl/2 

z     =     s   sin  P(l  -  hl  ; 

w     =     s   cos  0  +  T 

then  R  and  B  are  polar  coordinates  in  the  w-z  plane  and  the  mapping  of  the  shock 
(T  =  0)  and  the  wall  of  the  wedge  (s  =  0)  into  the  w-z  plane  is  given  in 
Figure  2.2.  The  shock  is  given  by  B  =  (3   where 


tan  p   =   (l  -  m^)1/2  tan  0 


and  the  wall  is  given  by  (3  =  0.   Since  0  and  iil  are  constants  the  w-z  plane  is 
a  linear  transformation  of  the  s-T  plane. 
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and 


Thus 


We  have 
3R 


3s 

£R 

3T 

9s 


£>T 


^  w 


Figure  2.2 


=     sin 


=     cos   P 


in  0(1  -  hl  )  '      sin  p  +  cos  0  cos  P 


1        c  pi  /p  "*) 

=     — =-    j  cos   P  sin  0(l  -  m_ )  '      -   cos  0  sin  P   £ 


R 


sin  p 


J^s,T;     " 


1        .      uf^  2xl/2 

-  —  sin  0(1  -  n^) 


so  the  transformation  is  non- singular  except  at  R  =  0  and  ©d  . 

We  will  look  for  solutions  that  have  a  singularity  at  R  =  0  which 
can  he  written  in  the  form 


I 


(51 


(R,  P)         Rntf^  (P) 
n=0 


^  =  1,    2,  ...,  6, 
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where  t  is  any  positive  constant,      jL  =m  ,       9 .    =  n,      £     =  \  and     £>  =  v. 


We  can  write 


V 


<5s 


■       E     R111'1    fntf  ^(sin  P  sin  0(l  -  m^)1'2  +  cos   p  cos  0) 
n=0  ( 

+  f    j(sin  0  cos  p(l  -  m..  )  '      -   cos  0  sin  p)  f 


P/. 


^     =       E       R1*"1    f  ntf^cos   0  -   fn/in  p} 
n=0  C  / 


*  df 

where         f     = 


d3 


We  assume  that  f   =  (m  )R_n  jf  0  and  substitute  the  above  relations  in  equations 
(1.2.11).  We  will  then  write  the  resulting  equations  as  power  series  in  R  and 
equate  each  coefficient  to  0.  On  substituting  into  ( 1.2,11)  we  have 

,nt-l 


E       Rn  "       <  ntf     j  cos   P  -  f      1  sin  p  "J 
n=0  (         n2/  n^  J 

- 

E   A^.R^"1  fntf  .(sin  0  sin  0(1  -  m?)1'2  +  cos  p  cos  0) 
n=0    J     C    nJ 

+  f  .(cos  p  sin  0(1  -  m  )  '   -  sin  p  cos  0)  7  . 

Now  A^.   is  a  function  of     r _,     £..     p       and    7v  and  so  can  be  written  in  the 

form 

00 
A   .      =       E       A     ,Rqt 
*  q=0        qVj 

where  A  ^  is  independent  of  R  and  is  a  function  of  the  fn^  f-,  ^  .  ••>  f  .  only. 

Therefore  we  can  write  the  flow  equations  above  as 
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Z       Rnt_1  ntf      [  Sp6  cos   p  -      Z     RqtA       ^.^(sin  p  sin  0(l  -  m2)1'2  +  cos   0  cos   p)] 
n=0  q=0 

=     Z     Rnt_1f      [(.sinp+     Z     RqtA      .<£,(sin  0  cos   Ptl-m2)1/2  -   cos  0  sin  p)] 
n=0  n        >e  q=0  qVJ    J6  ^ 

(2.2,2) 

where  the  summation  convention  is  used  for  Roman  suffices  running  through  1  and  2 
for  Greek  suffices  running  from  1  through  6. 

If  we  evaluate  the  term  in  R   we  get  a  set  of  six  total  differential 
equations  for  the  fy which  can  be  written  in  the  form 

6 


*ft  06 


where 

^V6  =  *v? 


^(L,  sin  P  +  k.tf'.Q.    (sin  0  cos  p(l  -  m  )  '   -  cos  0  sin  p)] 

» 

where   ayr  =  ^((^Jr^  +  ^±j^-r    <C^  and  ls  Siven  b^ 


~&ye 


Ell  -f05f03COSf04sin  P     °      ° 


f05(f^  '  1)cosf04Sln  P  "22 


E, 


E~  E  -sin  p     0       0       0 

\±  \2  0    -  sin  p    0       0 

1L_  I__  0       0      sin  p    0 

!?1  >2 

E6l  E62  0       0       0   -  sin  p 

—     —  2  l/2  — 

where  E   =  E  =  cos  P  sin  0(l-m  )  '  -  sin  p(cos  0  -  f  ^sinf  ,  ) 

and  the  other  E  ,'s  indicate  non-zero  terms  unimportant  to  the  discussion. 
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Now 
del 


et    |  |E   |  |   =  f  sin  0  cos  3(1  -  m^)1'2  -  sin  3(cos  0  -  fQI-   sin  f^H 


+  f05(l  "  f03}  cos2f04sin^ 

2 
>  0  since  f-.,  =  m      <  1  by  hypothesis. 

Since  in  addition  E   =  0  if  i  =  1  or  2  and  V=  3>  ^-,    5  or  6,  we  have 
f_  =  f  _^  =  0  for  all  3,  so  that  f-.  and  f   are  constants .   Since 

det||Ew|   =  sin  p  det  I  II.  .  I   >  0 

if  sin  3^0,  f  v  =   0  for  all  3  such  that  sin  3  =  0. 

Confining  our  attention  to  a  region  where  0  <  3  <  /7~we  have  at 
3  =  3n  (the  shock)  cos  0  =  sin  u  ,  \   =  f   =  1,  and  f  »  =  u   so  that 
cos  0  =  f nc-sin  f0K«  Making  this  substitution  into  (2.2.2)  and  rearranging  we 
get 

E     Rnt_1    [E  ,  f     -   ntB,    f       -   C    .]      =     0  (2.2.3) 

^  V&T16.  i/e  n<f         nVJ  \         ->/ 

n=0 


where 


/*  /  PI  /p 

^E^  sin  3  +  AQr  <b     (sin  0(1  -  n^)  '      cos   3  -   cos  0  sin  3)] 


.  pi  /p 

B>6  =    aw^4-<£  cos  ^  "  Aorj^^sin  ^  sin  0(x  ~  mi'       +  cos  0  cos  P)J 


and 


Cn2/   "  "a^r 


n       C'  2  l/2 

Z  A  _X.  V    (sin  0  cos  3(l  -  m,  )    -  cos  0  sin  3) 
|=1  q^]  J£[  n-q£v    r  L' 

+  (n  -  q)tf    As±n   P  sln  0(1  ~  mi)    +  cos  P  cos  0)) 

Thus  C  .1  is  a  polynomial  in  the  f~,,   f\^j  •••>  a^d-  f  -i^«  The  matrices  E  and  B 
are  given  below: 
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E  = 


P  1  /? 
(l  -  m  )  -  sin0cos3 


(-p  -  l)sin0sinP 

m 

1 

-[2  +  (/-l)m^]sin3 


0 


0 


sin  p 


-m  sin0sinp 


(l  -  m^)1^2sin0cosp      0       0 


0 


%k 


-  sinp      0 


0     -  sinp 


0 


0 


0 


0 


0 


E 


m. 


'25 


0 


0 


0 


0     -  sinp 


0 


0 


0     -  sinp 


P  1  /p      P 
where  E .  =  cos  0  sin  0  cos  3(1  -  m  )  '   +  sin  0  sin  3 

P  1  /p 
0  cos  p(l  -  m  )  '   +  sin  0  cos  0  sin  3. 


and    E   =  -sin 


B  = 


-  sin 


PI  /p      P 
3sin0(l  -  m  )  /c    -m  sin0cos3 


0 


( — -  -  l)sin0cos3 
-[2  +  (r-l)m^]cosp 


-sinPsin0(l  -  m  )  ' 


0 


0 


0 


0 


B 


2k 
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-  cosp      0 


0     -  cos3 


0 


m. 


cos  3 


0 


0 


cosp      0 


0     -  cosP 


where  B ,  =  -cos0sin0sinp(l  -  m  )  '   +  sin  0cos3 


.    2 


2  a/2 


and         B       =  sin  0sin3(l  -  m   )  /      +  sin0cos0cos3. 
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In  particular  Cnj/  =  C,^  =  0 

Now 

det||E^||  =  (1  -  m?)  sin20  sin  p  ^  0 

in  the  region  considered,  so  we  may  invert  it  to  get 

e 

f    =  nt  D  f   +  C  ,  (2,2,4) 

ny  V6    n  e    n  v 

from  each  term  of  (2.2.3 )»  where  the  C  ,  =  E  ^  C   are  also  polynomials  in 

-1 
sin  £c  cos  p,  the  f   and  f   for  m  =  1,  2,  ....  n-1,  and  the  matrix  D  =  E  B 
*  me.      me  7      '  '  ' 


and  is  given  by 


2 

ml 


,-     2J/2 


0        0        0        0 


[1   -  m1; 


2 


m2[2  +  m?(r-l)] 

cot6[2  +  (r-l)mj]    -± Xg   ,        cot  p     0 

(1  -  n±)  i 


0 


Vlk  P2k  0       cot  P     0 


D15  D  0        0       cot  P      0 


Cl      .  Ciml 


cot  p  i        0        0        0       cot  p 


(JUL  D — -  I- 

"l  (1  -  ^ 

(1  ~  m -J        2  ?  -,/p 

where   D  .   =   ^ tsin  0  +  sin0cos0cotp(l  -  m  )  '     ] 

"l 

P  p   n  /p 

D  ^     =     -[sin  0cotp  -   sin0cos0(l  -  iil)  '     ] 
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fi  2\!/2 

(1  -  nO        2  2  l/2 

D    =  -  p [sin  0cotp(l  -  dl,  )  /   -  sin0cos0] 

"l 
D    =  -[sin  0(1  -  m  )  '      +  sin0cos0cotp] 

In  particular  C  ,  =  C  ^  =  0 

The  solution  of  (2.2.4)  can  be  written  as  a  particular  solution  plus 
arbitrary  multiples  of  the  six  independent  solutions  of  the  homogeneous  equations 
obtained  by  setting  C  1  =  0„   Suppose  a  particular  solution  is  F  }   then  we 
may  write  the  general  solution  as 

f,  =  F.+a  F  ,  (2.2.5) 

nv     n>        n£  nVa-  v 

where  the  a   are  arbitrary  constants  and  the  matrix  F  ,,  is 
ne     ..  j  ny6. 


0 


2  2 

hl    sin  ntp  nLcosnt|3  0 


(l-m^)1/2cosntp        (l-m^)1/2sinntp  0 


F  „  F  „  (sinp)nt  0  0  0 


n31  n32 

'aJd-l  Fn^2 


*-m  F,,o  0  (sinp)nt  0 


Fn51  Fn52  0  0  (sinB)nt  0 

Fn6l  Fn62  0  0  0  (sinP)nt 
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2.3       Thus  we  have  expressed  the  ?.  as 

m  nt     2             2 

T,   =   E  R  (a  -.nLsin  nt0  +  a  „m,cos  ntp  +  F  .  ) 

^1       _  nl  1          n2  l          nl 
n=0 

>0     =       E    Rnt(-a  ,(1  -  mfh^cos  ntp  +  a  n(l  -  m^^sin  ntp  +  F  _) 
>2  „  v      nl  1  n2  1  n2 

n=0 

00 

E 

n=0 

In  this  and  subsequent  sections  we  will  discuss  the  boundary  conditions  for  a 
shock  attached  to  a  wedge.  This  means  that  if  the  flow  is  assumed  to  follow  the 


_       00  , 

J„  =   E  R^(a  _F  n  +  a  _F  /0  +  a  .(sin  p)nx  +  F  .),   2^=  3,  ^  5,  6. 
si/  _    v  nl  mJl  n2  nV2    nV x    rt  ny "       ->*   7  >-,» 


wedge  wall,  we  can  determine  a  ,  from  the  above  if  the  wall  is  |3  =  0.  The  next 
0      '  nl 

step  will  be  to  apply  the  Rankine-Hugoniot  equations,  which  can  be  written 

parametrically  as  f      =    <  (  £  )  on  the  shock  0  =  0  , 

This  gives  a  further  set  of  five  relations  to  determine  the  a  ,  and  t. 

n^ 

The  discussion  will  then  revolve  around  the  conditions  that  must  be  satisfied 
for  flows  of  this  class  to  exist. 

If  the  boundary  condition  at  the  "wall"  of  the  wedge  was  instead  that 
P  =  0  was  a  free  stream  line,  we  would  have  "?;(R,  0)  =  constant  so  that  a 
would  be  determined  instead  of  a     .  The  remainder  of  the  boundary  conditions 
would  be  unchanged.  We  will  only  consider  the  problem  where  the  flow  remains 
"attached"  to  the  wall,  but  it  should  be  noted  that  analogous  results  hold  in 
this  other  problem. 

We  will  discuss  the  problem  of  a  wedge  in  a  known  oncoming  flow  in 
the  R-3  plane.  The  shock  will  be  defined  as  T,  =  0  and  the  variable  s  will 
be  defined  along  it  such  that  X  =  1  on  the  shock  and  s  =  0  is  the  stream 
line  passing  through  the  apex  of  the  wedge  A.   In  this  section  we  will  assume 
that  the  flow  ahead  is  uniform  and  that  the  sides  of  the  wedge  are  straight 
for  a  finite  length.  A  shock  is  assumed  to  exist  that  will  turn  the  flow 
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Figure  2.3 

through  an  angle  V^  and  be  such  that  the  flow  behind  the  shock  is  subsonic  at 
A,  This  can  be  either  the  strong  shock  or,  when  the  angle  7^  is  large  enough, 
the  weak  shock.  The  boundary  conditions  that  must  be  enforced  at  the  wedge 
which  is  3  =  0  are  that  ¥       is  known,  or  in  this  case  that 


1, 


co 

=    z 

n=0 


*nv°> 


7T 


+ 


> 


or 


f^Co)  -  o 


n  >  1 


(2.3.1) 


If  we  chose  a  particular  solution  F  ^(3)  such  that   F  ^(0)  =  0  this  implies 
that  a       =  0   n  >  1.  This  is  the  only  condition  there  since  the  wall  can 
support  an  arbitrary  pressure  and  J    and  jp_  were  the  only  flow  variables  that 
were  necessarily  continuous  at  3=0.  At  3  =  3*  the  shock,  the  R-H  equations 
must  be  satisfied.  These  express  the  F      as  functions  of  a  =  JQ   -  rj.   Using 
a  Taylor  series  expansion  we  can  write 

j£,  =  fQ2,+  ay(a(s)  -  aQ)   +  0(a(B)  -  aQf  (2.3.2) 


30 

I are  functions  of  the  flow  ahead  of  the  shock  and  a^ 

only  given  by  the  R-H  equations  and  the  geometry. 
Now 

a(s)-a0    =   Jk-   J2  +  f0k-  f02 

=   Rt(f12O0)  -  f-^O0))  +  R2t(f22(P0)  -  f2i/e0)  +  ...        (2.3.3) 

Thus   since 

we  have  by  substituting  (2„3«3)  into  (2.3<>2)  and  equating  terms  in  R 

where  Q  is  a  polynomial  in  the  f  p(fO  -  f  ^(P/O  *"or  m  =  ®>  ^~*  2,...,  n-1. 
Since  we  know  ^ftj(P«)  le"t  us  assume  for  the  inductive  hypothesis  that  we  know 
the  f  i(Pn)>  m  =  0,  1,  . ..,  n-1,  so  that  Q.  is  known.  From  (2. 2'. 5) 

where  only  the  a   are  unknown,  so  we  have  on  substituting  into  (2.3o^-) 

or 


(2.3-5) 


if 


<VWV  +  FnW 

-  <W 

av 

"nsWV  +  Fnl«V 

-  V     " 

al 

2 
2hl    sin  2u 

ai  -      (r+i)      ^ 

0    . 

These  are  a  set  of  five  conditions  on  the  flow  for  each  power  of  R  ,  and  with 
the  boundary  conditions  (2.3.1)  at  the  wall,  are  the  only  conditions  that  apply 
in  the  neighborhood  of  the  apex  A  of  the  wedge. 
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We  will  also  need  the  following: 
Lemma   If  the  wall  of  the  wedge  is  straight  and  the  oncoming  flow  is  uniform 
and  if  f.,  (p  )  =  0  for  i  =  1,  2,  . ..,  j,   and  for  k  =  1,   2,      for  some  3  ,  then 

Xi£   X  X 

f.  ,(p)  a  0  for  i  =  1,  2,  ...,  j,   and  •  =  1,  2,  „..,  6. 

Proof   Assume  that  it  is  true  for  j  =  n-1.  Then  from  the  definition  of  C  / 
nv 

we  have  that  C  (  =  0  for  y=  1.  2.  ....  6     since  the  fn>,,  ....  f  , ,   are 
nV  '   '    '  16'    '      n-16 

identically  zero  and  hence  a  particular  solution  F  ^  is  identically  zero.  Thus 

from  (2.3-1)  and  (2.2.5)  we  have  a   _  =  0.  Then  since  f  _  (0  )  =  f  J$   )  =  0 

nl  nl  x     n2  x 

we  have  from  (2.2.5)  that  a       =  0.  The  Q,  of  (2. 3-5)  are  zero  when  the  lower 
order  coefficients  are  zero  as  it  is  a  polynomial  in  them,  so  we  now  have  the 
denominator  on  the  left  hand  side  of  (2„3o5)  zero.   Since  the  right  hand  side 
is  finite  we  must  have  the  numerator  of  the  left  hand  side  also  zero.  Thus  we 
have 

an^(sin  p0)nt  =  0      for  y=  3,  h,   5  or  6, 

or  a  1  =  0       since  sin  (3^0. 

Since  the  lemma  is  trivially  true  for  j  =  0  (there  is  nothing  to  prove  I)  we 
can  complete  the  proof  by  induction. 


Now  let  j  be  the  first  n  such  that  either  f  _(8)  or  f  A&)   is  non-zero, 

nl       n2 

From  the  lemma,  we  know  that  the  f  [$)   are  identically  zero  for  p  <  j   and  so  we 
can  write  from  (2.2.5) 


f ...   s=  (q:   sin  t.p  +  a_  cos  t.p)mn 
Jl       Jl      J     J2      jK/  1 


2 

(2.3-6) 

fj2  =   ("ajl  COS  V  +   aj2  Sln  t«5P)(1  "  4^      ' 
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From  (2.3-1)  we  will  have  a   =0  and  then  we  can  substitute  these  into 
(2.3.5)  to  get  either 

a2      ^   "  "V 


2        tan  tjP0 
1        m. 


(2.3.7) 


or  a.„     =     0 

J2 


Now  the  latter  is  not  true  for  otherwise  the  f  .,(3~)  and  f  _(p^)  would  both  be 

jl     0  j2     0 

zero  contrary  to  the  definition  of  j.     Thus    (2.3.7)  defines     t     to  be  one  of 
the  countable  set 

,n  2x1/2  >> 

t      .  1        cot'1!  ^  2  /    +     iF  (2-3-8' 

where  r  is  an  integer  such  that  t  >  0.  We  can  choose  any  such  r  and  any 

a_  and  in  general  get  a  solution. 
j2 

To  prove  this  we  must  show  that  the  remainder  of  the  boundary  conditions 

can  be  satisfied.  From  the  nature  of  the  solutions  (2.2.5)  the  conditions  (2.3.1) 

can  always  be  satisfied  by  choice  of  a  , .  For  v*  =  2   (2.3.5)  can  be  written  as 

nl 

2  a/2 


a2       (1  "  V    an2  Sln  tP0n  +  Rn2 


a  2 

1         nL.a  _.  cos  t  3^n  +  R  n 
J.  n2        0     nl 

where  R  n  and  R  _  are  known  functions  of  the  known  lower  order  terms  and  a  n . 
nl      n2  nl 

This  has  a  unique  solution  for  a  _  if 

n2 


fi    2xl/2 
ll     *   "  ml^ 


cot  t3Qn  /  -—   - =  cot  jtPQ   . 

2       m, 


Thus  it  has  a  unique  solution  unless 


BQt(n  -  j)  =  p7T 


for  some  integer  p. 
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This  will  be  satisfied  almost  nowhere  in  terms  of  the  oncoming  flow,  so  almost 
everywhere  there  will  be  a  unique  solution  for  a  -and  a  .  For  /^=  3>  ^$    5  or 
6  we  have  from  (2.3.5) 

where  P  ,  is  known  in  terms  of  a  . ,  a  _  and  the  lower  order  terms,  so  that 
nv  nl'   n2  ' 

since  sin  Pn  ^  0,  a  ,  is  determined  uniquely  from  these  lower  order  terms. 

The  flow  can  be  analytic  if  jt  is  some  integer  p.  Taub's 
necessary  condition  that  the  flow  is  analytic  is 

At,. a.   =  0       for  some  p 
2j  3 


or  that 


where 


ti  2\1/2 

a  (1  -  m  ) 

-~  tan  p9  ~ =  1   at  s  =  T  =  0     (2.3-9) 

2  m-. 


tan  9  =   (l  -  m  )    cot  u_ 


But 

2  l/2 
tan  (3   =  (l  -  m.  )  '   cot  u0  =  tan  9  , 

therefore  3  =  9  so  that  (2.3.9)  is  equivalent  to  (2.3.7)  with  tj  =  p. 
However,  this  condition  does  not  imply  that  the  flow  has  to  be  analytic  for  there 

are  many  choices  of  r  and  j  in  (2.3.8)  that  will  make  jt  other  than  an  integer. 

9  a 
The  curvature  of  the  shock  is  -= —  which  behaves  asymptotically  like 

s     as  s  — *  0  along  the  shock.  Therefore  the  curvature  will  be  finite  but 

non-zero  at  A  only  if  tj  =  1.   If  tj  <  1  then  the  curvature  will  be  infinite 

at  A,  if  tj  >  1  it  will  be  zero. 

When  (2.3.9)  is  satisfied  with  p  =  1  the  shock  angle  is  known  as 

[k] 
the  Crocco    angle.  There  we  have  from  (2.3.8) 


tj    =    1  + 


3^ 
r7T 


^0 


Since  Pn  <  7T/2  and  tj  >  0,  we  must  have  r  >  0,  so  that  the  curvature  at  A 
will  be  finite  if  r  =  0  and  zero  if  r  >  0. 

For  shock  angles  between  the  Crocco  angle  and  the  extreme  angle  we 
have  tj  =  a  +  rT/$       where  0  <  a  <  1.  Thus  if  r  =  0  the  shock  has  infinite 
curvature  at  A  but  if  r  >  0  it  will  have  zero  curvature  there. 

For  shock  angles  between  the  Crocco  angle  and  the  sonic  angle  we  will 
have  tj  =  a  +  rT/p  where  1  <  a  <  TJ2.,  Since  (3n  <77/2  and  r  is  an  integer 
we  must  have  r  >  0.  Therefore  tj  >  1  and  the  shock  has  zero  curvature  at  A. 

We  have  shown  in  this  section  that  in  general  a  countable  choice  of 
r  and  j  in  (2.3.8)  and  an  arbitrary  choice  of  a   will  lead  to  a  set  of  the  J 
which  satisfy  the  boundary  conditions  at  the  shock  and  locally  at  the  apex  of 
a  straight  wedge.  These  parameters  must  in  some  way  be  determined  by  the 
boundary  conditions  along  the  sonic  line  that  starts  from  the  shoulder  of  the 

wedge. 

[11] 
We  also  rederived  the  necessary  condition  of  Taub    for  an  analytic 

flow  to  exist  and  saw  that  this  was  not  a  sufficient  condition. 

2.4       Now  let  us  consider  the  problem  of  a  wedge  in  a  known  non- uniform 

oncoming  flow. 


7-7  7  7  /  /  /  /  /  /  7 
Figure  2.k 
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Let  AC  "be  the  wedge  side  and  XA  be  the  known  stream  line  passing  through  A.  The 
Rankine-Hugoniot  equations  will  give  a_  in  terms  of  ~/s  and  the  oncoming  flow. 
This  is  assumed  to  exist  and  be  such  that  the  flow  behind  the  shock  at  A  is  sub- 
sonic. We  assume  that  the  flow  ahead  of  the  shock  is  analytic  everywhere  con- 
sidered, that  the  shock  is  analytic  except  at  A,  and  that  the  flow  behind  the 
shock  is  analytic  except  at  A.   The  suffix  0  refers  to  the  flow  ahead  of  the 
shock .  Now 

\  -  (JA'-r^ti.^c-2-!)] 
9h    (X 


so 


hrr  do- 


9s 


dS/  °  Y+  1  +  2rt<r    -   1)        9? 


Now 

U  v  y 

r=   -^    =    c^^p-V       and        uto    =    Vln(i2o-ao) 

therefore 

ds  cQ-   cQ     ^  9s JO  v0cQ      Us        '        2      Qs 


'0 


hence 


%. 


=  ka_  +  B,  (2.4.1) 

9s       1    1 

where  B  is  a  function  of  the  flow  variables  ahead  of  the  shock, 

only,  and  so  is  dependent  on  the  shock  direction  but  not  its  derivatives.   If 
(D  is  the  rotation  ahead  of  the  shock,  and  L  the  specific  enthalpy  we  have 
(see  Taub1   equations  (4.17b)  and  (4.19a)  ) 
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nO 
so  that 


„      i    =v^+c2 


O^s  ~0 


^s/<y      n0,s     "     Un0^   "   C0^s/0     +  ^-1 


2co    ^0 
2s 


B, 


k<r 


3s  /0+  ^+l  +  2«<r2-l)  c0    Ws/°      v?= 


W<%)      +UnO 


Unoi"-co(3i)o 


0  0 


-^^-^[v-^-ao(Do] 


(2.4.2) 


Now 


J2     =     k+£,     =  1 


=     ka2  +  B2 


u     -J&  -  u+  V 

n      3s  t       ps 


(2.4.3) 


where  B     is  also  a  function  of  the  flow  ahead  and  the  shock  direction  only. 
We  look  for  a  solution  of  the  form  (2.2.5)  behind  the  shock.     At 
3=0,   the  wedge,   we  have 

co 


jf  (T,0)     =     0     =       E       T^f      (0)    . 

n=0 

Let  us  assume  that  the  shape  of  the  wedge  is  a  known  analytic  function  of  T. 

Then  we  must  have     t  =  l/p  for  some  integer     p  >  0  and     f      (0)     =     0,   unless 

p    |   n   .      Since     1  <  g-    <  m     ,      a     >  p  "    ,      >  0.      Thus  we  get  from  (2.1.4), 

(2.4.1)  and   (2.4.3) 


3s    "  Bi 


ka. 


a. 


ka. 


at 


T  =  0       if  k    4     0 


JT'B2 
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or  that 


z  Rn/p-l  _5_    (   )(1  .  ^Bln2   )V2  .  Bi 


P 

Z  R       T"  fn2(P0)(l  "  V111  °!0)    "  B2 


'2 


Taking  the  limit  as  s  or  R — »0  and  assuming  B  ,   B  /  0 

lf  JL.    >     i  ,  (a) 


a2       B2  p 


a^     B±   -  f  .1QQ)(1  -  n^sin2^)1/2      ^  _^ 
a2      B2  -  f  2(Po)(l  -  n^sln2^ 


if  -J-  =  1      (b)    (2.^.10 


or 


a 


~   =    J     ta\  if  -sL  <   1  (c) 

a2      fJ2(PQ)  P 

where  j  is  the  first  n  such  that  either  f  ,(6^)  or  f  ~(p«)  is  non-zero.  Now  if 

nlv  0     n2  0 

j  <  p  then  we  have  from  the  lemma  of  section  2„3>   f .  (p)  =  0  for  i  <  j,  so 


that 


F~  p  +  aj2  ' 


f ._  =  m_  (a.,  sin  -*—   p  +  a.^  cos 
,ll      xv  ,il     P 


f '   =  (1  -  m?)1'2^.,  cos  -^-  p  +  a.,  sin  -&-   )   . 
j2         TL'      jl      p  K    j2      P 


But 


f.o(0)  =  0,     so  a.n  =  0, 
,l2v        '  ,il    ' 


so  from  (2.^.4) 


2 


a£      (1  -  ^) 


Thus  if  there  exists  a  positive  rational  number  r  <  1  such  that 


cot(rP-)  = 


38 
ft        2.1/2 

1     (1  -  m£) 


0'      a_        2 
2       ^ 

there  can  exist  a  solution  of  this  type.  However  this  is  not  in  general  true 
so  we  must  have  j  >  p.  Then  we  can  group  cases  (a)  and  (b)  of  (2.^.^4-)  together 
to  get 


al      V1  "  *^*ta£%Y        -   fpi^o^ 


l2      B2(l  -  m^sin2a0)'1/2  -  f  (pQ) 


We  will  have 


f  .,  =  m,  (a  nsin  £  +  a  ^cos  p) 
pi     ±N  pi        p2 

f    =  (1  -  m^n^-a  cos  3  +  a  sin  p) 


and. 


„  /.    2v-l/2 
apl  =  -KgU-n^) 


where 


3=0, R=0 


2 


so  a  _  is  the  solution  of 
P2 


t>   ft  2.2      v-l/2  2  .  __     2    .      .    (-  2v-l/2 

ll  l'      "   misin  a<y  "  ap2micos   ^o  +  ^i3111  ^      ~   "V 


a«  -a    ft  2    .    2      x-1/2  /n  2*1/2    . 

2  B2(l  -  Yin  <V  "      p2^      "   *V        Sin  P0  "   ^C°S   P0 


E2- 

P2( 

(2.4.6) 


This  exists  and  is  unique  if 


2 

a  m 

2       (1  -  m£) 

«.  •   .  [^  1 

This  is  violated  at  the  Crocco    angle,  so  that  at  the  wedge  angle  which,  in 

the  case  of  a  straight  wedge,  allows  an  analytic  solution  we  have  no  solution 


39 


to  this  more  general  problem  unless  (2.4.6)  is  satisfied  when  a       =  0. 
(2.4.7)  is  violated  again  when  ~y^  =   0  which  implies  that  a     =    7l/2  =   3  and 


that  a  =  0.  Then  we  must  have 


B1  =  -K2m1 


(2.4.8) 


If,  in  addition,  flow  ahead  is  irrotational  and  is  isenthalpic,  then  we  have 
from  (2.4.2) 


B, 


'Jh\ 


?s 


4m, 


0     y+ 1  +  2  ri^Q  - 1) 


U 


X. 


1  /  ^1 


n0       2c, 


<?S    /0 


Unq!o  /  % 
2      f    .3 


or 


P 


B. 


I. 


-m 


0  V  DT 


2m£  -    (r+  3) 

0     2^  -  (r  -  1) 


Then  (2.4.8)  is  precisely  the  condition  of  Tsien    .   If  (2.4.8)  is  not 
satisfied  then  a  solution  of  this  class  does  not  exist. 

2.5 

to  an  inequality  when  a  shock  involving  a  singularity  of  the  type  S~'  "'  is  used. 

We  have  seen  that  it  is  not  one  of  the  family  of  solutions  discussed  and 


ro  ] 
Lin  and  Rubinov    discuss  a  solution  which  relaxes  Rsien's  condition 

,1/2  , 


I  believe  that  such  a  solution  does  not  exist,  so  let  us  discuss  their  technique 

[3] 
for  finding  a  solution.   Similar  remarks  apply  to  the  work  of  Chow    in  his 

extension  of  the  results  to  the  psuedo- stationary  problem.  For  the  normal 

shock  problem  let  us  suppose  we  know 

/  -  u  =  a  =  -£-   +  bs^  +  o(s^)   for  <£>  0 

where  b  and  6   are  parameters  to  be  determined.  We  can  then  calculate 


3} 

Tr        (3,0)   =   Bl- 


ktffs2*'1 


r+ 1  +  2r(m£  - 1) 


/  2c(-1n 
+  o(s    ) 


'Jt 


<?s 


(s,0)     =     B2   + 


S&>/~\±  -   m?) 


1+0 


2  +  (r  -  iK 


+  o(s       )  . 


Now  we  can  write 


c? 


?Ji 


-/i(s'T)  =  A(s>0)  +  T    ?k  (s,0)  +  0(T } 

=  ^.(s,0)   +  TA±J     -y^   (s,0)  +  O(T^) 


where 


(A..)     .     =      (1  -  m2)"1 
ij   s=0  v  1 


0 


1  - 


nu 


in  this  case.  Therefore 


/2(0,T)   = 


JL  +  T 

2      2 
ml 


k/^2 


m. 


>mQb 


lim  s 


2*- 


^+  1  +  2<T(ii1q  -  l)   s-?0 


1  -  M  +  0(T2) 


Since 


(2.5.1) 


^2(0,T)  =  -f-  +  iy?  +  O(T^) 


.2^-1  .. 


we  must  have  either  6  =   l/2  or  B  =  -K^m   so  that  lim  s     is  finite 

s^O 
and  non-zero.  Thus 


B, 


2m' 


S2 


\ 


m. 


m^    r+   1  +  2^(niQ  -  1) 


or  0 


(2.5.2) 


which  gives  the  Lin  and  Rubinov  condition  that 

Km2  +  B   >  0  . 

Since  a  shock  that  behaves  like  s  '   around  s  =  0  requires  that  j/p  =  l/2  in 
the  treatment  of  section  2.k   we  must  satisfy  equation  (2.^.5)  if  "the  Lin  and 
Rubinov  solution  is  to  be  a  member  of  the  family  of  solutions  discussed.   However 


since  a  =  0  and  P  =  TJ2}    (2.4.5)  is  not  satisfied.   Therefore  I  do  not 
believe  that  a  solution  of  the  type  they  suggest  exists.  The  false  reasoning 
in  the  above  outline  of  the  argument  must  be  that  the  error  term  0(T  )  of 

equation  (2.5.1)  is  not  bounded.   This  can  happen  since  it  was  the  limit  as  s 

-1/2 
tended  to  zero  of  an  earlier  term  -which  may  behave  like  s    ,  and  so  become 

infinite . 

Now  the  error  term  0(T  )  can  be  expressed  as 

-~  T   — =~  (s,0)        for  some  Q6.  (0,T)   . 

We  cannot  evaluate  this  explicitly,  nor  get  a  lower  bound  for  it,  but  we  can 
show  it  tends  to  oo  at  Q  =   0  as  s  tends  to  0. 

p/  l 
+  0(s  °   )  at  T  =  0 


*7e 

1            ^>-2 

2 
2^ 

2             2 
1  -   m        3  s 

1      ^"VCT  -  1)  2D(1  -  ng) 


1  -  m^      (T-  1)  iHq  +  2 


+  Ofs^1) 


which   — »  od   as  s  — >  0  for  o  =  l/2.   Consequently  we  have  no  guarantee  that 
(2.5.1)  is  a  representation  of  J     at  s  =  0. 

If  we  estimate  the  radius  of  convergence  of  J      expressed  as  a  power 
series  in  T  about  the  point  (s,0)  we  will  have 


00 


n   ^ 


72(s,t)  .  E  -E-  -^(x,o) 

n=0  n!   P  T 
The  radius  of  convergence  is 

r  1  •       I  1          ^2  |l/ni-l 
L  lim  sup   I  — - —  I   J 

n-^oo      n°   pTn 


k2 


For  n  even  we  have 


where 


Sa- 


°(^n) 


,n     2j   3sii     n2 


/"n2b(l  -  m?) 


>  =   ^(n)  .  01     +  o ;(/-n} 


Ss 


n 


(r-  1)  iEq  +  2 


where  <((n)  =  /(/-  l)  --•■(/-  n  +  l). 

^n-<      ^(n)n,   cJ-n/-.    2x 


Therefore 


n/2 


c9T 


n 


(r-  l)m£  +  2 


m ,  -  ±  I 


*   o(s^n) 


Therefore  the  radius  of  convergence  may  be 


< 


lim  sup 
n-»co 


4^  i""  i.i  * 


\  1/2 


1  -  m 


1 


-1 


=   (1  -  m1)    |sj 

if  d     is  not  an  integer.   Thus  the  region  of  convergence  of  these  power  series 
may  be  as  given  in  Figure  2.5  in  the  neighborhood  of  A. 


Shock 


Figure  2.5 


^3 
2,6      A  note  on  pseudo- stationary  flow. 

The  equations  of  psuedo- stationary  flow  can  be  written  in  the  form 
(see  Taub    ) 

•J  %  d} 

9}  3}  2}        r 

2     y   1  'r   ^1  *Z         ?    2     5     y 

5r  -  h{  y3  - 1]  cos  h  tt  +  -vlnn  tt  -  y^;  C0SA 

Tr   •   -fe+(r-i>J3l   7*"    jj 

%       x  sh 

Tt    =    cosA   55" 

Tt    --j^  Tt    "    x  (2-6-1) 

or  if      |  | a. . | |      can  be   inverted  as 


*J„  ^ 


?F    =    Vji+Hv  <2-6-2' 


where 


4l2  ^  COS  JM      "A22^  COS  JV    "  ^  "   V2  +    (^-   ^jV'    V    °'    -1  "   Hl  f3 

We  will  make  the  transformation  to  the  R-f3  plane  and  substitute 

J„     =     Z    Rntfn  O) 

n=0 
as  before.   H  .which  is  a  function  of  the  J     can  be  written  as 


hk 


oo 


Hy         HO)  RDt 

q=0   H 


where  the  H  ,  are  functions  of  the  f  ,  for  m 
qV  mc 

differential  equations  will  be 


=  0,  1,  2,    ...j    q.   The  resulting 


E  Rtn_1(f  i  -  1U  f  nt  -  C  ,)  =  Z  B^H  (ft)  Rqt 

v  nv    Vc  ne-      nv         >6  q£v 
n  q 


(2.6.3) 


In  order  to  match  terms  we  must  require  that  for  each  q  >  0  there  exists  an  n 
s  ich  that  nt-1  =  qt.  Thus  t  is  the  reciprocal  of  an  integer  p.   If  we  substitute 
for  this  t  into  (2.6.3)  and  equate  terms  in  R     we  will  get 


nV  p   >6   n£    nv 


(2.6.k) 


where 


Jn*  "  CnV 


=  C  ,  +  E~Jh 

nV    ^6  n-p6- 


0  <  n  <  p 
p  <  n 


In  particular 


V*W 


m  sin  ft 


ft  2\        ' 

_         c    (1  -  m   ; 


-    cos    ft 

:i  -  m?)1/2  -J 


1c1(l  -  n^)" 


(2.6.5) 


so  that   if  C     ,  is   zero,    a  particular  solution  of   (2.6.^)  for     n  =  p     is 


V(p)  = 


o, 


sin  ft 


-X. «     JC  •     X.  o     a 


Yi(1  ~  mi) 

We  will  prove  that   if  a  solution  for  r  of 

/i  2xl/2 

ll  '     "  ml' 


(2.6.6) 


cot  rft 


a. 


a. 


2 


(2.6.7) 


^5 
is  a  positive  rational  number,  then  a  solution  to  the  local  boundary  conditions 
in  regular  reflection  can  be  written  in  the  form 

CD      / 

r,  =  T,      Rn/P  f  .0)  where  r  =  -J-  . 

Jy  rt       nv  p 

n=0  * 

In  particular  the  first  n  such  that  all  of  the  f   (p)  are  not  identically  zero 


is  n  =  min(j,p). 

First  note  that  Taub     has  shown  that  an  analytic  solution  exists 
if  (2.6.7)  admits  an  integer  solution.   Since  the  class  of  functions  we  are 
considering  includes  the  analytic  functions,  it  is  already  proved  for  integral 
r,  so  we  only  have  to  consider  those  r  which  can  be  written  as  j/p  where  p  /[  j 
(p  does  not  divide  j).  The  problem  of  regular  reflection  is  equivalent  to  the 


Uniform 

Pseudo- stationary 

Flow 


/////////////// 

Figure  2.6 


problem  of  a  wedge  in  a  uniform  pseudo- stationary  oncoming  flow,  so  we  will 
consider  this  latter  problem  only.  The  flow  behind  the  shock  at  the  apex  of 
the  wedge  is  assumed  subsonic. 

The  boundary  conditions  at  the  wall  behind  the  shock  are  that 
f  ?(0)  =  0  for  n  >  1.   At  the  shock  we  have  (see  Taub     equations  8.2  and 


8.8) 


*t±  952     cos  ?i 


-     kan , 


5s  l*        as       yG 


and  we  can  express   the  other  conditions   similarly  as 

Sgj    =     k(s)  av(s)  +bv(s)  (2.6.8) 

where  a,(s)  and  b^(s)  are  functions  of  the  >  (s)  but  not  their  derivatives. 
Nowk(s)R  can  be  expressed  as 

ClRt  +  C2R2t  +  °°°° 
and  both  a^(s)  and  b^(s)  can  be  expressed  in  the  form 

Dt      „2t 

ay  +  a^   R  +  a  „R   +  . . . 

where  the  C  ,  a  ,  ,  and  b  ,   are  functions  of  the  f. ,   i  =  0.  1,  ....  m  only. 

m   >m      Pm  16       *   *    '      " 

Further  we  note  that  the  a ,  ,  C  and  b  ,  are  zero  if  p  /  m  and  if  all  the 

vm.       m      ^m  ' 

f .(B)  are  identically  zero  for  all  1  <  i  <  m  such  that  p  /[  i.  This  happens 

because  the  a  ,  R  ,  C  R   and  the  b  ,  R  '  are  sums  of  products  of  the  f  R 

rm.       '     m  >m  nfe 

with  coefficients  which  depend  on  the  f_.  only,  so  that  they  can  be  written  as 

sums  of  terms  like  77"  f   ,  (8)  where  E  m.  =  m.   Each  term  in  the  sum  is  thus 

m.  e.  .i 

ill  i 

zero  unless  p  divides  each  m. .  However  this  implies  that  p  divides  m  contrary 
to  hypothesis. 

Thus  if  we  take  the  terms  in  R     from  (2.6.8)  we  will  get  equations 


of  the  form 


f  .    a  ,_C  +  R  ; 
n^      *>n    nV  (2^9) 


f       an_C  +  R 
nl      10  n    nl 


where 


n-1 

R  .  =   Z   C  a_,     +  bv 

nV  q  V  n-q    y  n-p 

q=l   ^ 

and  thus  are  functions  of  the  lower  order  f  ,.     Therefore  the  R  are  also  zero 

m£  n> 

if  p  /  n  and  if  all  of  the  f  ,  are  identically  zero  for  1  <  m  <  n-1  such  that 
1  m6  —   — 


P  /f  m. 


hi 

Returning  to  the  differential  equations  (2.7.^-)  we  note  that  by  the 

same  reasoning  the  J  ,  are  identically  zero  under  the  same  conditions  that  the 

R  ,  are  zero. 
nV 

Now  we  are  in  a  position  to  set  up  the  inductive  hypothesis.   Suppose 

all  of  the  f   are  identically  zero  for  m  <  n  -  1  such  that  p  j  in  and  p  y[  n. 

Then  if  n  <  j,  the  f  .  are  identically  zero. 

Note  that  the  R  ,  and  J  ,  are  both  zero  under  the  conditions  of  the 
nv      nV 

hypothesis  so  that  a  particular  solution  of  the  differential  equations  (2.6.^) 

is  zero.   In  the  general  solution  (2.2.5)  we  have  that  a  ,  =  0  from  the  con- 

nl 

dition  at  the  wall  of  the  wedge.   Then  from  (2.6.9)  we  have 

a2       (1  -  mj  a  sin  n  PQ/p 

=  o s~~ 7 —  (2.6.10; 

a_  2a  Ocos  n  pn/p 

1  m.  n2       0' 

so  that  if  n/  j.  a  ^  =  0.  Then  the  denominator  of  the  left  hand  side  of 
*   n2 

(2.6.9)  is  zero,  so  that  since  the  right  hand  side  is  well-behaved,  the  numerator 
of  the  left  hand  side  will  also  be  zero.   On  substituting  the  general  solution 
from  (2.2.5)  we  will  get  a     ■  =  0  for  y=  3,  k,    5  and  6.   Thus  we  have  shown 
that  the  a     ,  are  identically  zero. 

Thus  by  induction  we  can  show  that  the  f  ^(@)  are  identically  zero  for 
all  n  <  j   such  that  p  /[  n.  When  n  =  j   (2.6.10)  admits  an  arbitrary  solution 
for  a   since  r  =  j/p  satisfies  (2.6.7).  For  n  >  j  there  will  be  extra  terms 
in  the  equations  (2.6.10)  so  these  will  then  have  a  non-zero  solution  for  the  a. 
in  general. 

The  cases  that  are  still  unsolved  are  thus  those  where  there  does  not 

exist  a  positive  rational  solution  of  (2.6.7).   A  method  for  finding  a  solution 

in  these  cases  is  to  express  7  as  Z  E  R    f   (S).   Substituting  into  (2.6.2) 

-J*  nmv 

n  m 

we  get 


m 

f   ,  -   (n  +  mt)f   Dv.  +  C  i 
nmv    v        nm£  y*         nm* 

where  the  matrix  D  is  as  before  and  the  C  j  are  polynomials  in  the  f    and 

nmy     r  *  rsifr 

• 

the  f  t   for  r  =  0,  1,  . . • >   n  and  for  s  =  0,  1,  . . . ,  m,  except  that  they  do 
rs  ^ 

• 

not  include  the  f   or  f  .   In  particular 
mn     mn 

W  ■  ° 

V  ■  ° 

and  C   is  the  E^  Hn^  given  by  (2.6.5).   By  techniques  similar  to  the  treatment 
of  stationary  flow  we  can  write  the  solutions  as 

1  =  f001  +  M^Ol8111  P  +  ai02  COS  ^  R  +   ml^a011Sln  tP  +  a012  COS  tP^R  +  "• 

J2  =  f002  +  (1  "  4)l/2("ai01COS  P  +  ai02Sln  P)R 

+  (1  .  4)V2(     ^  tp  +      in  tp)Rt Rsin  ,  +  _ 

m1c1(l  -  n^) 

(2.6.11) 
Since  -— = =  =  0  at  (3  =  0  we  have  a    =  a    =  0.  From  the  R-H  equations  we  have 

^(R.Pq)  =  f001  +  a1(a(s)  -  aQ)  +  0(a(s)  -  aQf 
and 

y  COS   (J.  / 

J2(R,P0)  =  fQ02  +  a2(a(s)  -  aQ)  -  -—±   (l  -  m2  cos2  j^)'1'2!*  +   ... 

These  can  be  obtained  directly  by  integrating  equations  (2.6.8)  and  the  fact  that 

o   p   1  /p 
R  =  s(l  -  m^cos  u  )  '        at  £  =  £  .   Substituting  these  into  equations  (2.6.11) 

and  using  the  fact  that 


1*9 

a(s)  -  aQ    -     (f012(P0)  -  t01k%))^  *  (floe(P0>  "  WP0»  R  +  ••• 

we  get  on  equating  terms   in  R     and  R 

-m^012cost30     =     a2(f012OQ)  -  f0ll+O0))     =     (1  -  4)1/2a012sin  tBQ 


al 


(2.6.12) 


and 


mia!02COS   P0     =     a2(f102(P0)  "   WP0)}  =   ^  "  ml>         ai02COS   P0 


al 


(2.6.13) 


since 


o    n  In 

cos  p.  cot  u.  tan  B   (l-m   )  sin  P 

77        2  '     2      vI72     =     "7.         7.  2l     ~2      7T/2     =  Sec   8  =     T-,      2*1/2 

(1  -  m  cos  p.,  )  '  (1  +   (1  -  m   )cot  u    )  '  sec  Pq  (l-"1-,  ) 

i      -if*!    d-4)1/2? 

Since  we  are  assuming  that  — r—  cot  s ^ 7      is  not  an  integer  we 

P0      (  a2      m\  J 

have  a-,np  =  0  from  (2.6.13).  t  is  given  by  (2.6.12)  and  o:m   is  arbitrary. 

Thus  once  again  we  have  a  coutable  choice  of  t  and  an  arbitrary  para- 
meter to  choose,  both  of  which  must  be  determined  in  some  way  by  additional 
boundary  conditions .  The  remaining  parameters  will  be  determined  uniquely  in 
terms  of  these  provided  that 

a         (1  -  m   ) 

cot[(mt  +  n)B0]     /     -±-    ~- =     cot  tp 

2       hl 

for  all  n,  m  >  0  except  for  n  =  0  and  m  =  1.  Again  we  can  say  that  this  happens 
for  almost  all  oncoming  Mach  numbers  and  wedge  angles  since  it  requires  that 

B0(n  -  l)t  +  m  /  p7T  (2.6.110 

for  some  integer  p. 
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/i   2\l/2' 
1    ^.i(a!   d-^) 


2.7      We  can  raise  the  question  of  what  happens  if  — —  cot   < ^~ 


is  a  rational  number  and  we  use  a  double  series?  A  more  general  question  is  what 
would  happen  if  we  used  a  multiple  series  of  the  form 

n,t-  +  . .  .in.  t 

J    =      E  E  ...  E  R   x    ^  k  tn,...^0)  (2.7.D 


'v 


nl  n2    nk 


in  the  stationary  case  and  of  the  form 

n,tn+. .  ,+n.  t,  +m 

?  =  E  ...  E  E  r   x    ^k   f      (p.)    (2.7.2) 

nl    nk  m  '    k. 

in  the  pseudo- stationary  case?  If  we  follow  the  techniques  of  section  2.6  we 
will  get  a  number  of  equations  of  the  form  of  (2.6.12)  and  (2.6.13)  each  giving 
rise  to  an  arbitrary  parameter  a  if  the  t.  is  a  solution  of 

0      O      O    l—  1 


cot  t.P 


/,  2xl/2 

ll        (1  "  ml} 


10  aQ  m2 

d  m 

Thus  we  can  have  a  solution  of  this  form  with  k  arbitrary  parameters  if  the 
higher  order  boundary  conditions  can  be  satisfied.   They  can  if  there  do  not 
exist  positive  integers  n.,  at  least  two  of  them  being  non-zero,  such  that 

-rr 
n,^  +  .  ..+  n^  =  t±   +  -g2-  (2.7.3) 

in  the  stationary  case.   This  will  not  be  satisfied  in  general  since  we  can 

IT 

write  t.  =  t,  h r —  q.  where  q.  is  an  integer,  so  that  (2.8.3)  is  the  con- 

i    1    p  n 

dition  that 

(nx  +  n2  +  ...  +  n^t^  -  t^  =   (p  -  ^VV^n^     (2.7.*0 

77/(3  t   is  not  a  rational  number  in  general  so  that  this  does  not  admit  a  solution 
with  n  +  n  +  . . .  +  il   >  2.   In  the  pseudo-stationary  case  the  equivalent 
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condition  is  that 

(nx  +  n2  ...  +  nk)t1PQ  -  t^  +  m^Q     =   mpQ  =  %  -  ^...n^)   (2.7-5) 

which  is  not  satisfied  in  general. 

Thus  there  are  equivalent  solutions  to  the  stationary  and  pseudo- 
stationary  problems  of  the  form  (2.7.1)  and  (2.7.2)  respectively,  where  the  t. 
satisfy 


cot  t.p 


/i    2^/2 
i±        (1  -  m1) 


I  o      a_,       2 

2       m 

There  are  k  arbitrary  parameters  provided  that  equations  (2.7.^-)  and  (2.7«5) 

cannot  be  satisfied  in  the  respective  cases. 

In  conclusion  it  is  worth  noting  that  there  are  numerous  techniques  that 

could  be  applied  to  this  and  related  problems.   For  example  in  the  preceeding  a 

particular  class  of  orthogonal  functions  of  the  independent  variable  R  was  used, 

namely  a  set  of  different  powers  of  R.   If  )  0  (r);  n  =  0,  1,  ...  {    are  a  set  of 

orthogonal  functions  of  R  such  that 

d0 

-  (R)  =   F  (0„(R),  ...  L    (H)) 


dR   v  '  "   n^0x  "     rn 

and 

0  (R)0  (R)  =  G  (0n(R),  ...  0  ^  (R)) 

where  the  F  and  G  are  linear  functions  of  the  0's,  then  the  use  of  a  series 

^        n      n 

Z  0  (R)f  <(3)  for  the  f.  would  result  in  a  set  of  linear  total  differential 
n=0 
equations  in  the  f  ^((3)  whose  coefficients  are  functions  of  the  f nj(&)>    ••• 

f  -lyCP)  only.  We  could  ask  what  set  of  orthogonal  polynomials  lead  to  the 

simplest  solution,  or  if  any  other  set  will  yield  a  solution  to  some  of  the 

remaining  unsolved  problems.  More  generally  we  could  ask  are  there  suitable 
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independent  variables  with  which  to  use  this  technique  in  partial  differential 
equations  of  the  type 

ik    -  A  „  *2± 

for  arbitrary  matrices  A  ? 
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III  VISCOUS  FLOW 

3.1      The  algebraic  R-H  equations  and  the  algebraic  matching  conditions 
across  non- viscous  shocks  and  slip  streams  only  admit  trivial  solutions  when 
applied  around  the  triple  point  in  Mach  reflection  and  no  solution  when  applied 
to  regular  reflection  if  the  incident  shock  is  weak  and  the  angle  of  the  inci- 
dent shock  close  to  the  angle  at  which  the  flow  behind  the  shock  is  sonic. 
Experimentally,  however,  both  regular  and  Mach  reflection  can  be  observed  in 
this  region  (see  Smith    )„   Since  the  techniques  of  Chapter  II  required  that 
shocks  exist  that  turn  the  flow  through  the  appropriate  angles,  those  techniques 
cannot  be  used  in  discussions  of  Mach  reflection  and  regular  reflection  for 
which  no  solution  to  the  R-H  equations  or  no  non-trivial  solutions  to  these 
equations  exist. 

The  incident  shock  being  weak  implies  that  the  ratio,  of  the  pressure 
on  either  side  of  the  shock  is  close  to  1  or  that 

Pl      .    2T   ,2   .  v  _  . 
=  l  + -  (a-    -1)^1 

p0  y+ 1 

2 
so  that  o  ££l.   Then  if  the  flow  behind  the  shock  is  nearly  sonic  we  have 

from  (1.2.1)  -  (1,2.4) 

2       (ag  -  c-2)(r+  D2<r2  +  ((r-  D<r2  +  2)2 

m. 


1         (2  +  (r-  l)cr2)(2r<r2  +  1  -   V) 

2 
Since  a-   ^>1  we  get 

1  -  ml  ^  mo 

so  that  the  cases  of  greatest  interest  are  those  for  which  the  Mach  number  of 
the  oncoming  flow  is  close  to  1. 


We  are  going  to  try  to  modify  the  R-H  equations  and  the  matching 
conditions  across  stream  lines  by  introducing  viscosity  and  then  considering 
the  result  'as  viscosity  tends  to  zero.  The  equations  of  viscous  flow  are 
given  by  (l.l.l)  through  (1.1.3).   We  assume  that  the  viscosity  w  is  a  constant 
If  we  replace  the  line  integrals  in  (l.l.l)  -  (1.1.3)  by  area  integrals  of  the 
divergence  we  can  then  express  the  equations  as  second  order  partial  differential 
equations  in  four  unknowns.   These  are  too  complicated  to  be  solved  analytically 
and  until  more  is  known  of  the  appropriate  boundary  conditions  for  a  given 
problem  it  is  difficult  to  see  how  to  integrate  them  directly  numerically. 

We  can  transform  the  independent  variables  x  and  x  of  equations 
(l.l.l)  -  (1.1.3)  into  the  s-T  coordinates  defined  in  section  1.2  to  get 
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=  0 


(3.1.1) 


"1   rS. 


T0  S0 


(  (OU  u.  +  \.<f.  .p  +  wH.  )ds   - 
M  n  1    j  l/T     1 
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r   1    1 
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t 

J       (p 

•J 
V 

£..    +  wG. )dT 

U          1      J 

=    0 
(3.1.2) 

and 


"1  rs 


1  2 

^n  -i- 


TQ  *0 


1    rTi 

+  yfj  pUn  +  vl)dsl  +  [^  J     vJdTJ  =  0   (3.1.3) 
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H. 


2Yij        ,  ^v  #2  ^2    v 
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In  problems  involving  the  matching  of  conditions  across  discontinuities  in  non- 
viscous  flow  one  uses  the  conservation  equations  around  two  types  of  contours, 
first  around  contours  like  PQRS,  where  PQRS  contains  a  shock  discontinuity 
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Figure  3.1 
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which  is  a  discontinuity  across  which  matter  flows,  and  secondly  around  contours 
like  ABCD,  where  ABCD  contains  a  stream  line  discontinuity  across  which  no  matter 
flows.   In  each  case  one  lets  the  region  shrink  to  the  line  of  discontinuity, 
and  then  by  appealing  to  the  fact  that  the  result  is  true  for  arbitrary  lengths 
of  these  discontinuities,  or  by  differentiating  with  respect  to  the  length, 
arrives  at  the  conclusion  that  the  jump  in  certain  quantities  is  zero. 

To  simplify  the  viscous  problem  we  assume  that  there  will  be  regions 
where  viscosity  is  important  corresponding  to  the  shocks  and  slip  streams  of 
non- viscous  flow,  but  that  elsewhere  the  viscous  terms  will  be  negligible.   Thus 
we  will  consider  the  flow  pattern  shown  in  Figure  3^2  as  the  viscous  analogy  of 
the  non- viscous  pattern  of  Figure  3-1-  We  can  shrink  contours  like  PQRS  and 
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Figure  3-2 
ABCD  until  PQ,  RS,  AB,  and  CD  are  the  bounding  lines  of  the  adjacent  non- 
viscous  regions,  and  so  have  relations  between  the  flows  at  PQ  and  RS  and 
between  the  flows  at  AB  and  CD.   KLMN  will  be  the  intersection  region  where 
two  or  more  transition  regions  interact. 

In  the  non-viscous  problem  the  definition  of  the  various  lines  like 
PS,  QR,  DA,  BC,  etc.,  presented  no  problem  since  as  the  limit  was  taken  they 
shrank  to  zero  and  the  result  was  independent  of  their  choice.   In  the  viscous 
case  the  complexity  of  the  results  may  depend  strongly  on  their  choice. 
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The  contour  PQRS  is  analogous  to  a  contour  containing  a  shock  dis- 
continuity, so  we  will  call  a  viscous  transition  region  a  shock  region  if 
matter  flows  in  one  boundary  between  the  region  and  the  adjacent  non-viscous 
region  and  out  of  the  other.   Thus  PQ  and  RS  are  boundaries  between  the  viscous 
and  non- viscous  regions.  We  will  define  PS  and  QR  as  stream  lines  since  these 
are  independent  of  the  choice  of  coordinates. 

We  will  classify  a  region  as  a  slip  stream  region  if  no  matter  flows 
across  the  boundaries  between  it  and  the  adjacent  non- viscous  regions.   Thus  if 
ABCD  is  a  slip  stream  region,  AB  and  DC  will  be  stream  lines,  or  lines  of 
constant  s.   To  simplify  the  treatment  we  will  define  AD  and  BC  as  lines  of 
constant  T.   Since  we  can  expect  a  region  of  this  type  to  originate  at  some 
disturbance  such  as  the  intersection  region  KLMN  we  can  define  T  =  0  to  be  one 
boundary  of  KLMN,  for  example,  KN.   To  define  KN  relative  to  the  flow  we  must 
look  for  characteristics  that  distinguish  the  intersection  region  from  the  slip 
stream  region.   In  this  case  we  may  expect  the  intersection  region  to  have  some 
of  the  characteristics  of  a  shock  region,  namely,  large  derivatives  of  flow 
variables  along  the  stream  lines.  We  do  not  expect  these  in  a  slip  stream 

region,  so  a  definition  of  KN  is  the  boundary  between  regions  where  the  viscous 

£  v 
forces  due  to  terms  such  as  w  .=-=  are  significant  or  are  not  so. 

<y   1 

One  further  simplifying  assumption  we  will  make  is  that  the  "shock 
thickness"  is  a  constant  along  its  length,  or  that  the  length  of  the  stream 
lines  within  the  shock  region  is  constant.   This  means  that  if  T  =  0   is  defined 
to  be  the  boundary  RSNK,  then  QPML  is  also  a  line  of  constant  T. 

The  effect  of  the  assumptions  will  be  to  make  viscous  effects  due  to 

terms  in  w  —  dominant  in  a  shock  region,  while  viscous  effects  due  to  terms 

o'T  °   ' 

9v 

in  w  —  will  dominate  in  a  slip  stream  region.   The  intersection  region  will 

have  characteristics  of  the  adjacent  shock  and  slip  stream  regions. 
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We  have  tacitly  assumed  continuity  conditions  across  "boundaries  which 
we  should  explicitly  enumerate.   Along  stream  lines  we  will  assume  that  all  flow 
variables  are  continuous.  Across  stream  lines  we  must  have  r     continuous  as 
this  is  a  geometric  variable.   Since  we  are  taking  account  of  viscous  forces 
which  arise  from  the  derivatives  of  the  velocity  we  must  also  consider  the 
velocity  to  be  continuous  and  differentiable  (the  derivates  need  not  be  con- 
tinuous) so  that  across  stream  lines  the  velocity  is  continuous.   If  the 
derivatives  of  velocity  have  jumps  across  stream  lines,  the  viscous  forces  will 
have  jumps  so  that  we  may  expect  the  pressure  to  have  discontinuities  also. 
The  pressure  discontinuities  will  arise  at  the  boundaries  of  stream  line  regions 
because  of  the  approximations  to  the  flow.   If  we  could  solve  the  equations 
exactly  they  would  not  arise. 
3.2      When  we  examine  a  shock  region,  supposed  bounded  by  the  two  lines 

T  =  T„  and  T  =  T  ,      we  differentiate  equations  (3.1.1)  -  (3.1.3)  with  respect 
to  s  in  a  similar  manner  to  the  method  of  non-viscous  theory.   The  resulting 
equations  are 
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The  lines  T  =  T   and  T   can  be  considered  as  belonging  to  the  non- viscous 


regions,  so  that  the  viscous  terms  on  the  left  hand  sides  of  (3.1.1)  -  (3.1.3) 
are  zero;  therefore,  the  remainder  of  the  left  hand  sides  are  the  set  of  R-H 
equations  equated  to  terms  of  the  right  hand  sides  instead  of  to  zero.  Thus 
the  right  hand  sides  represent  correction  terms  which  must  be  estimated  to  the 
accuracy  required. 

To  estimate  these  terms  we  must  have  some  knowledge  of  the  flow  in 
the  interior  of  the  shock  region.   In  order  to  understand  the  structure  of  the 
region  we  will  examine  a  shock  region  where  all  functions  are  independent  of  s, 
that  is,  the  region  is  one- dimensional.   We  will  study  a  linear  approximation 
to  the  flow  in  this  region  and  then,  in  Section  3-3;  extend  this  approximation 
to  a  shock  region  where  derivatives  with  respect  to  s  are  small  but  not  zero. 

If  the  shock  is  one- dimensional  the  right  hand  sides  vanish  so  we  are 
left  with  the  R-H  equations  alone.   To  the  approximation  that  the  flow  is  non- 
viscous  outside  of  the  shock  region  these  are  correct,  but  we  have  not  yet 
determined  the  thickness  of  that  region.   To  do  that  we  must  also  consider 
equations  (3.2.1)  -  (3.2.3)  evaluated  at  various  interior  points  of  the  shock 
region.   Then  we  get  the  equations 

(°V    =  M  =  f>0vQ  (3.2.4) 

Mv  +  p  +  —  w  -~-  =  MvQ  +  pQ  (3.2.5) 

2  v2 

MV^7^I+  "T  wv^  "  "-§-♦*<*>  T^T  (3.2.6) 

which  have  the  solution 

T  =   3M^+lKv^-1)    [V0  1O^0  "  V)  "  Vl  l0g(v  "  Vl)]      (3<2*T) 

f   _  JL 
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(vQ  -  v)(v  -  v1) 


p  =  MvQ  -  pQ  -  Mv  +  M(f+   1)  — .  (3.2.9) 

The  flow  variables  p.  P  and  v  attain  their  extreme  values  p  ,  fi     and  v  or 
p  ,  p     and  v  at  T  =  -  oo  or  +  oo  respectively.   To  simplify  the  solution  we 

could  use  (3.2.7)  -  (3.2.9)  for  part  of  the  plane,  say  from  T  =  p-  to 

T  =  — —  and  a  non- viscous  flow  in  the  remainder  of  the  plane.   This  would  corres- 


pond to  picking  points  where  the  flow  had  attained  say  957"  of  its  final  value,  or 
any  other  suitable  value.   Another  common  definition  of  shock  thickness  is  to 
relate  it  to  the  total  jump  in  a  flow  variable  divided  by  the  maximum  value  of 
the  variable's  derivative  in  the  T  direction.  This  corresponds  to  taking  that 

variable  to  be  linear  and  fitting  it  at  the  point  of  maximum  slope.   Both  of  these 

aw 
definitions  lead  to  a  shock  thickness  of  the  form  — for  weak  shocks,  where 

mQ  -  1 

m  is  the  oncoming  Mach  number  and  a  is  a  constant.   This  requires  that  a  solution 

of  the  flow  equations  be  known,  as  indeed  one  is  for  the  one-dimensional  shock. 

With  the  two-dimensional  flow  in  mind  we  look  for  approximations  to  the  flow  in 

the  one- dimensional  problem  that  can  be  defined  in  terms  of  the  flow  equations 

(3.2.4)  -  (3.2.6)  only.   The  simplest  approximation  we  can  make  which  takes 

account  of  viscosity  is  a  linear  one  for  p  and  v  where  for  T  lying  between  Tn  and 

T  =  T 0  +A    we  write 

T  -  T 

p  =  P0  +  (P2  "  p0}  ~7T- 
v  "  vo  +  (vi  "  V  "J— 

M 


f    = 

V 

These  satisfy  (3.2.4)  so  we  must  choose  A    to  try  to  satisfy  (3.2.5)  and  (3.2.6) 
in  some  manner.   Since  in  this  approximation  (3.2.5)  is  always  in  error  by 
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—  (v  -v  )  we  can  choose  A    to  satisfy  only  (3.2.6)  at  one  point,  say  the  mid- 
point T  =  T +A/2.      If  a  higher  order  approximation  were  to  be  used  we  could 
satisfy  both  equations  at  more  points.   Substituting  T  =  T  +A/2     into  the 
approximations  and  evaluating  (3.2.6)  we  get 

i6w(v0  +  V]L)(r  -  1) 

^  =   3(  V  -  1)  M(vQ  -  v±) 
From  (1.2.1)  -  (1.2.4) 

V^  -  Vn 
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/  w   \ 

so  the  thickness  of  the  shock  in  this  approximation  is  0  /  — for  weak 

2  (  m0  "  V 

shocks  where  m  %il.  > 

3.3       In  this  section  we  will  discuss  the  nature  of  an  approximation  to  the 
structure  of  a  two-dimensional  shock  region  similar  to  that  used  for  the  one- 
dimensional  shock  region  in  the  previous  section.  The  approximations  used  inside 
the  shock  region  are 


\ 


=   a 


+  bv(s  -  s.)  +  cv(T  -  T_  +  dy(s  -  s.(T  -  T  J     y=  1,  2,    ...,5 


where  X    ,   c  .  ^  ,  ^5",  and  \   are  p,  y  ,   v,  \  and  a  respectively.  Equation 
(3.2.1)  specifies  jo    in  terms  of  these  other  variables.   T  =  T   and  T  =  T 
are  the  two  boundaries  of  the  shock  region.  These  forms  for  the  variables  will 
be  assumed  for  T  between  T_  and  T  only.  We  have  already  assumed  that  the  thick- 
ness of  the  shock  is  a  constant.   Let  us  assume  that  the  equations  have  already 

aw 
been  used  to  determine  it  in  the  form  — . 

mQ  -  1 
We  will  substitute  the  relations  (3«3.l)  into  the  conservation  equa- 
tions (3.2.2)  and  (3.2.3).   In  addition  we  will  substitute  them  into  the  two 
geometric  equations 
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J-  f^  ^J 

\~]     =   I     l  cos   |i     — |     dT  (3.3.2) 

L      _         /  Q  s 

rp  J  rn 

0  0 

[*}=]  ^    fs     ^  .  (3.3.3) 

m  m 

xo  0 

We  are  assuming  that  the  flow  ahead  of  the  shock  region  is  known  and 
that  the  "shock  angle"  a  is  known  on  one  boundary.   We  will  take  X   =  1  on  this 
boundary.   These  conditions  specify  ten  of  the  twenty  unknowns  of  (3.3. l). 

The  five  equations  (3.2.2),  (3.2.3),  (3-3.2)  and  (3.3-3)  are  then  func- 
tions of  s  which  would  vanish  identically  if  the  approximations  (3«3»l)  were 
exact  solutions.  We  can  pick  the  remaining  ten  unknowns  of  (3.3«l)  to  satisfy 
the  five  equations  at  s  =  s  and  the  derivatives  of  these  five   equations  with 
respect  to  s  at  s  =  s  .   Thus  the  flow  behind  the  shock  region  is  specified 
in  terms  of  the  flow  ahead  of  the  region,  the  angle  of  the  shock  region  boundary 
and  its  first  derivative  with  respect  to  s.   If  we  had  taken  a  quadratic  in  s 
instead  of  a  linear  approximation,  the  flow  behind  the  shock  region  would  also 
have  depended  on  the  second  derivative  of  the  shock  angle. 

We  are  going  to  examine  these  equations  in  some  detail  when  the  flow 
ahead  of  the  shock  region  is  nearly  sonic  since  this  will  simplify  the  equations 
and  this  is  the  case  in  which  we  are  most  interested  as  was  mentioned  in  Section 
3.1. 

We  will  assume  that  the  Mach  number  of  the  flow  ahead  of  the  shock 

region  is  given  by 

2     ,1/2 
m0  =  1  +  —  6 
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and  that  the  angle  of  the  boundary  T  =  T  of  the  region  is 

T 
a     =  -~-  +  af 

where  a  is  of  order  1  in  £   but  may  depend  on  dr.     We  will  denote  differentiation 
at  the  value  of  s  under  consideration  with  respect  to  a  by  ' .and  differentiation 
with  respect  to  s  by   .  Thus  we  have 

k  =  &    =     b.6 

Let  us  take  the  flow  ahead  to  be  uniform  with  A.  =  p_  =  1.   The  R-H  equations 
have  the  solutions 

3  "  £_         f-i    _      2\a2\    J.  csl^ 


\     =    ^(1  -  2(r/l)   (1  "  &)6)   +   °(f  } 


P  -  1  f^;  (1  -  a2)^2  +  o(e^)  (3.3-10 

Ut  =  -,/fat(l+-^-  (3  -  a2)^2)  +  0(e   5) 
and  the  corresponding  derivatives  are 

P  =  -  —73  a^k  +  0(<fT) 

\     =  &   JT^L     a^k  +  0(^7)  (3.3.5) 

2 
U+   =  -/Tk  -  -f-  /?  (1  "  a2)k  +  0(6  8) 


2+ 

where  k  has  been  assumed  to  have  order  £    .  The  reason  for  this  assumption  is 

that  it  can  be  shown  to  be  consistent  with  later  work  on  two  and  three  shock 
region  intersections.   In  addition  we  will  assume  that  the  shock  thickness  A 
is  of  order  6        (this  is  in  agreement  with  the  results  of  one- dimensional  flow) 
and  that  the  solutions  of  the  equations  for  the  linear  approximations  are  of 
the  form 


6k 


u  =  u  -  -^—  (Tte*  +  Vt*   )  +  o(*  °) 

n      n    V  +   1  / 


/ 
f 


P  +— fi  (fe2  +V06k)   +  0(6  6) 


ut  =  ut  +^3  +  0(^5) 

(  (3.3.6) 


u 

n 


U   =  u  +  v,'k62  +  0(6  8) 

where  if/,   0,  ^  . . .  etc  „  are  independent  of  6-   . 

When  we  evaluate  the  Integrals  below  we  can  use  the  trapezoidal  rule,,  as 
errors  will  have  orders  in  £  at  least  as  high  as  those  to  be  neglected  in  other 
terms  evaluated.    It  may  be  verified  that  when  equations  (3.2.2)  and  (3.2.3) 
are  used,  the  viscous  terms  wG.  and  wJ  will  have  a  higher  order  in  £.   than  any 

a    fTl 

terms  that  are  kept.   It  will  be  verified  for  the  term  dr-  \   wJdT  in 

cps  J 

0 

Appendix  I,  the  proof  for  other  terms  is  similar. 

We  will  need  -—— •  evaluated  behind  the  shock  region. 

IF  1   u+ 

J2     =     a  +   u  =  -g-  +  a£  +  tan    -fi- 
ll 

Substituting  from  (3.3.6)  and  from  (3.3.^)  we  get 

*r-  ^     (  a6  +  4-  (3  -  a2)A  -  -hL  &3        ) 
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Differentiating  with  respect  to  s  we  get 


3* 


£.   k     ~>    ,  akd- 


2 


8, 


■yf      ^=r  -  — ^~r     k^2(l  -  3a2)  ->"-:-  -V  '  ^—  +  0(^u) 


(3.3.7) 
We  must  evaluate  \   and  a  ahead  of  the  shock  region.  From  (3.3.2)  and  (3-3.3) 
we  have 


V 


;i  »x 


6 


>x 


2    gs 


+  o(^6) 


(3.3.8) 


and 


sin(  J  -  a)dT 


=  0(65) 


(3.3-9) 


We  apply  the  conservation  equations  (3.2.2)  which  can  be  written  as 


>Uu.  +  p\.<*  .    =  -~ 


o> 


pu 


1J 


£ 


dT  + 
J    v         1 


"0 


£ 


(3.3.10) 


where  the  6  .    are  the  viscous  terms  of  orders  &      and  6      for  i  =  1  and  2 

i 


respectively.  The  proof  of  a  similar  result  for  the  energy  equation  is  dis- 
cussed in  Appendix  I„  Multiplying  (3.3»10)  hy  \         and  summing  on  i  we  get 

1 


U, 


c\ 


=    \ 


A    -A   (pu,  %   L  +  o(d5), 


iT    2  ~   <ps  x*  J   v   "T 


"G 
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on  using  the  trapezoidal  rule  to  evaluate  the  integral  and  using  the  fact  that 
the  flow  ahead  of  the  shock  is  independent  of  s. 
The  right  hand  side  can  be  written  as 

iTx   2   v   v  *  £s    j   v   £s 

=  4(psin^-  cosn.^  )  +  0(tf  5) 


J?s 


3s 


since  \   =  1.   From  (3.3-7)  -~  =  0(.r- 7)  and  from  (3.3.*0  H  =77/2  +  0(6   ). 

Therefore  this  can  he  written  as  -  -^ ^  +  0{  6    ).   From  (3.2.1)  the  left 

hand  side  is 


Wf. 


-1    ut- 


\ 


1       U, 


i-  (pVt0  [  —J  +  Wt0  ^3 


/r?< 


by  (3-3.^).   Therefore 


7' 


2/F  ^S 


(3.3.11) 


Differentiating  with  respect  to  s  we  get   V  =  0  since  — £  =  q  hy  the 
linearity  assumption.   Multiplying  (3.3.10)  by  d1.  .\.T  we  get 


U 


?\ 


n 


\ 


+  p\ 


"0 


=  ^ikW  -4-  £■  (pu.  -^  )  +  ou6) 

1        2      ^  s    vxr  j      v 


or 


Tj 


"0 


u2 


f 


+  P 


p(l  -   \) 


+   ^ik^kT. 


A_ 

2 


-^-   (pu.     -M)   +  0(66) 
^)s    VJr   J         V 


By  use  of  (3 • 3 .8)  the  right  hand  side  can  be  written  as 
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% 


is 


Vz) 


f-     IT    T,  +  (  sln  "  -gi  +  p  cos  v-  —sj\       2 


A  +0U6) 


aci     -4"  +  0(<f6)      • 


On  substituting  from  ( 3 . 3 « 6 )  and  using  (3.2.1)  we  get 


r 


fi(>+  V?*^2+7T1  (0^062)^2+-|--fI  f{±-  *2)±h 


r 


a^r 


&  -£   ^  +  0(66)   . 


(?s    2 


(3.3.12) 


The  energy  equation  (3.2.3)  may  be  written  as 


v  n 

M  -j—  +  ^p 


U 
] 

r+ 1 


=  o(*°) 


"0 


6 


since  the  viscous  terms  have  orders  £     and  higher,  as  is  verified  in  Appendix  I, 

2  2    2    2. 
Using  X   v  =  U  +  U7  and  (3-3.8)  we  may  rewrite  this  as 
n     u 


1       u2-^2 

m  -a s  x_ r 


u 


2\ 


/-    1   *        *  J 


Y 


f1pUn(  JL-   1)]   +  0{6G) 


_  _  JWT    J^    %  +  o(6e) 

-     r-  i      2     js  +  ul'    j    • 


Substituting  from   (3.3-6)  we  obtain 
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f  -  1  \  r    J        2<r+D  (r-D2(r-i) 

(   (mf(f-D  (r-ifif-D     2 

-  7^7  4  ^t  +  0(e6)    •  (3-3-13) 

From  both  (3-3.12)  and  (3-3. 13)  we  find  by  equating  terms  in  6:     that  0  =  \J-\ 

This  is  valid  only  at  the  value  of  s  we  are  considering.   It  is  not  a  functional 

relation  because  we  are  unable  to  satisfy  the  equations  identically  in  s  with 

this  approximation.   In  addition  we  insist  that  the  derivatives  with  respect  to  s 

of  (3.3.12)  and  (3-3.13)  are  satisfied  for  this  value  of  s .   It  can  be  seen 

immediately  that  they  will  lead  to  relations  0'  =  HJ    at  this  value  of  s. 

Multiplying  (3-3.13)  hy  — — -   .  subtracting  (3.3.12)  and  substituting  for  V| 

Y  n  [ 

from  (3.3.11)  we  get 

if  g\   *&*   ££    y.(l-  a2)  k         A       *£  (**ik) 

Substituting  for  ~ —  from  (3 • 3 • 7 )  and  for  %^-     from  (3-3.6)  we  get 

(4a  -  0')La  =  -|-  (L  +  y)(f   +  2(1  -  a2))+  0(62)        (3.3-15) 

where  L  =  6r      \A  .      We  must  also  equate  the  differential  coefficients  of  (3.3-12) 
and  (3.3.13)  to  0  at  the  value  of  s  under  consideration  to  get  two  further 
relations  between  1^ ,  y>' ,   0  and  0'.  We  have  already  remarked  that  this  will 

lead  to  0'  =  T6-' .   Since  (3.3.1U)  is  a  linear  combination  of  (3. 3.12)  and 

52/      32 
(3.3.13)  we  can  differentiate  this.   Noting  that =  ^  =  0  by  the  linear 

3s  ds 

approximation  we  get 

-  y-0'  -  (1  -  a2)0'  +  2a^  +  -|-  (ka  -   0' )   =  0(dr2) 
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or 


(4a  -  0')(l  -  a2  +  V-+-|>)  =   (2(1  -  a2)  +  /)2a  +  0(^2).     (3-3.16) 
One  solution  of  (3.3-15)  and  (3.3.16)  to  order  1  in  6-   is 

0'   =  4a,     7*  =  -2(1  -  a2)   . 

This  corresponds  to  p  =  1,   U  =  U  .,  or  in  other  words  there  is  no  change 
through  the  shock  region.  The  non-trivial  solutions  are 

f  '  '-f-^^t  A^"+  -ty2+kLa2}  ■        (3-3-i7) 

o 

If  k  =  0  then  L  =  0.   The  solutions  are  then  ~f  =   0  or  -(l  -  a  ) 

and  01  =  0  or  co  respectively.   Since  the  latter  solution  would  involve  infinite 
derivatives  of  velocity  outside  of  the  viscous  region,  it  will  be  discarded. 
When  L  is  non-zero  it  is  not  obvious  which  branch  of  the  solution  should  be 
used.   We  will  use  the  branch  corresponding  to  the  zero  solution  at  L  =  0. 
The  flow  variables  can  now  be  expressed  by  equations  (3. 3-17)  and 

0'        =     -—^    [/+  2(1  -    a2)][/    +  L]   +  Ji-a  (3-3.18) 

p     =    i  +  yn  [2(1  "  a2)  +  t]&2  +  °(&k)  (3.3.19) 

2       =    ^(1  +  f^i   C2(l  -   a2)  +  >  ]  ^2)  +  OU  k)  (3.3.20) 


c 


v2       =     r(l  -  -yn   [2(1  "  &2)  +    f]  &2)  +  °(  ^  (3.3.21) 

-  2     =   —  "  —  a(1  "  a  }  *    "  7TT  "  2(r+i)  ( !a  ~  0  ]  +  0(  ~  ) 

(3.3.22) 

f*        =     yf-i     k^(0'    -    4a)   +   0(^T)  (3.3.23) 


<2J 


2 


f  +  1 


k(l  -  3a2)^2 


r 
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(3.3.2)+) 


These  will  be  referred  to  as  the  modified  R-H  equations.   The  correc- 
tion terms  depend  only  on  L  which  is  proportional  to  kw.   We  will  assume  that  it 
is  sufficient  to  introduce  these  extra  terms  only,  and  study  the  effects  of  these 
modifications  on  problems  which  have  no  solutions  when  the  unmodified  R-H  equa- 
tions are  used. 

3^ 

A  B 


Viscous 
Flow 


D 


Figure  3.3 

We  will  next  study  the  viscous  analogue  of  a  slip  stream.   In  non- 
viscous  flow  we  applied  the  four  conservation  equations  to  a  region  ABCD 
containing  a  discontinuity.  When  the  region  was  shrunk  to  the  line  of  dis- 
continuity we  obtained  equations  that  had  a  solution  where  the  line  of 
discontinuity  was  stream  line,  so  that  there  was  no  mass  flow  across  the  dis- 
continuity,  and  that   [pu.  — "■]  =  0  across  the  stream  line.   This  implied  that 
J"    and  j    were  continuous  across  the  stream  line. 

In  viscous  flow  we  define  the  region  between  AB  and  DC  (Figure  3«3) 
as  a  slip  stream  region  if  the  boundaries  AB  and  DC  are  stream  lines  and  if 
the  flow  above  AB  and  below  DC  is  non- viscous  to  the  approximation  used.   Across 
AB  and  DC  we  assume  that  ^   and  v  are  continuous.   Continuity  of  X  is  not 
assumed  as  was  stated  in  Section  3.1°  p  is  defined  to  satisfy  the  conservation 
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of  mass.   We  apply  the  conservation  of  momentum  and  energy,  equations  (3.1.2) 
and  ( 3.1. 3 )>  where   s  =  sn  is  CD  and  s  =  s..   is  AB  and  then  differentiate 
with  respect  to  T  to  get 

1    &  -      rsi 

pu.  -ii    =  -g-       I    (/jUu.  +  p\.<f.  .  +  wH.  )ds  (3.^.l) 


J   v  J 


S0  S0 


and 


where  the  left  hand  sides  are  evaluated  in  the  non- viscous  region  so  that  the 
viscous  terms  are  zero  on  these  sides „ 

Using  the  linear  approximations  (3.3°l)  for  each  variable  we  get  a 
set  of  relations  involving  the  coefficients  of  (3.3.1)  which  could  he  satisfied 
identically  in  T  if  the  approximations  were  exact  solutions.  We  will  enforce 
the  equations  at  a  given  T  =  T   and  also  ask  that  the  T- derivatives  of  (3°J+»l) 
vanish  at  T  =  T..  .   If  the  viscous  terms  are  small  the  integrand  in  (3.^.2)  is 
approximately  the  enthalpy,  a  quantity  which  is  nearly  constant  along  a  stream 
line.   We  want  to  enforce  this  over  as  many  stream  lines  in  the  viscous  region 
as  possible.   With  a  linear  approximation  we  can  only  enforce  it  on  two  stream 
lines  which  we  will  take  as  s  =  s   and  s  =  s  ,  so  we  have 

—  vV  Un  -| —  +   •  1  pUn  +  wl)  =  0     at  s  --  s0  and  s±        (3.^.3) 


Sl     ■  .  , 
Equations  (3.  .l)  and  its  T  derivatives  express      pu  „  — - 

Sl       U.  -,  so 


and 
L-  "  j   v  _j 

•  1  ° 

p  ~:  .  .  — "—  I   in  terms  of  the  flow  in  the  interior  of  the  slip  stream 

region.   If  the  flow  along  AD  is  known  then  we  can  use  (1,2.9)  and  ( 1.2.10)  to 
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find  — =   and  — —  at  A  and  D.   (3. ^-.3)  expresses  — £  inside  the  slip  stream 
region  in  terms  of  the  other  variables.  The  continuity  conditions  give  ""7  and 
v  inside  the  slip  stream  identical  to  J     and  v  in  the  non- viscous  regions  at  the 
boundaries  AB  and  DC,  so  that  (3A.I)  will  give  a  set  of  four  relations  between 
the  flow  either  side  of  the  slip  stream  region  and  p,  X   and  a   along  AD.  We  are 
going  to  examine  these  relations  in  some  detail  for  a  "weak"  slip  stream  region, 
which  is  a  slip  stream  region  with  the  following  properties :   If  the  flow 
variables  are  scaled  so  that  at  some  point  in  the  flow  m  =  p  =  o  =1  and 
J>p  =  "/2  then  throughout  the  region  we  have 

c2       v2  2 

p,  f   ,  m  —  and  —  have  the  form  1  +  0(  6     ) 

?2  =   -f  +  0(^3) 
a    =    JL   +  o(^) 

X     =     1  +  0(  6     ) 

or  in  other  words  they  have  the  same  magnitude  as  the  variables  behind  a  weak 
shock  region.   In  addition  we  will  assume  that  the  T-derivatives  of  the  variables 
along  AB  and  DC  in  the  non- viscous  region  have  the  same  magnitude  as  the  T- 
derivatives  behind  a  weak  shock  region,  namely 


)  >        -jr  =  u^  er'  *     > 

-3 


I  =  0(6k)   ,    3I  =  O(^)      , 


that  the  thickness  of  the  slip  stream  region  D  =  s   -  s   =  0(^      )  and  that 

h 

the  pressure  discontinuities  at   s  =  sn  and  s  =  s   are  of  order  £    .   These 

assumptions  will  be  shown  to  be  consistent  for  regular  and  Mach  reflection 
later.   We  define  A,  =  1  at  T  =  T,  or  AD,  and  since  a  is  a  linear  function 
on  AB,  we  define 
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a   =  aQ  +  (a^  -  a0) 


s  -  s. 


on  AD. 


Let 


a. 

1 


IT 


b.6  +  0(  e5)  i  =  0,  1. 

i  ; 


Substituting  into  (3oU.l)  we  get 


pu 


ij 


J  v 


S0 


"5T 


!l   /  n   Ui  ^   ,n       ^2     Ui  ^2   _   )p  >, 


r  5     ij. 

where  the  o  .  are  the  viscous  terms  and  have  orders  £   and  6  for  i  =  1  and  2 

i 

respectively.  A  proof  similar  to  that  in  Appendix  I  will  verify  this.   Also  we 
can  verify  that  the  viscous  terms  in  (3«^-3)  can  be  neglected.   Substituting  for 

-=r-=  from  (3-^-.3)j  neglecting  terms  of  order  o  .  and  using  -«r-=  =  -  ,-.  m 

from  (1.2.11)  on  the  non- viscous  boundary  we  get  for  i  =  1 


\ 


f{  bl(%)Sl-o(^Js0^ 


E 


* 


^)8l^^-)-\#)^^-)-^ 


oj 


+  0(£5) 


(3.^A) 


and  for  i  =  2 
■l 

P 


'-^1^ +  o(£U) 


^1 


where  I  -^-=  I    is  evaluated  in  the  non-viscous  region  on  AB 


&i 


~75~r  i  "*"s   evalua"t'ed  in  the  non-viscous  region  on  DC, 


7^ 


Differentiating  (3.^.1)  with  respect  to  T  and  substituting  we  get 

Sl_ 


pu. 


J   v 


S0 


=  € 


where  the  6 .  represent  the  terms  on  the  right  hand  sides  of  the  T-derivatives 
of  (3.^.l)  and  have  orders  £•   and  £  respectively.  Thus  we  have 
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=    o(£7) 


(3.M) 


and 


"5~t 


=  0(6°) 


(3.1*. 7) 


for  i  =  1  and  2  respectively.   Equations  (3«^-«^)  -  (3«^-«7)  are  four  relations 


between  p  and  5   on  each  boundary  of  the  slip  stream  region  and  the  flow  along 
T  =  T  or  AD.   The  only  parameters  on  this  boundary  are  b  ,  b  and  D,  since  v, 
L,  p,  \  and  p  are  known  as  linear  approximations  from  the  continuity  con- 
ditions and  the  additional  assumptions.   These  assumptions  were  that  the 
variables  had  certain  orders  which  will  be  shown  to  be  consistent.   If  the  stream 
lines  are  straight  and  the  pressure  is  constant  along  them,  then  the  right  hand 
sides  of  (3.^.^)  and  (3 •^••5)  are  zero.  The  left  hand  sides  then  give 


soL  J 


=  0  , 


which  are  the  matching  conditions  across  a  non- viscous  stream  line. 
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3.5       We  have  discussed  linear  approximations  to  shock  regions  and  slip 
stream  regions.  Now  we  will  consider  two  problems  in  which  shock  regions  inter- 
act and  produce  a  slip  stream  region.   Such  interactions  occur  in  regular  and 
Mach  reflection  of  shocks  shown  in  Figure  3«^.  They  also  occur  in  other 
instances.   In  Figures  3-4  the  region  ARMB  will  be  referred  to  as  the  inter- 
section region.  The  line  RM  will  correspond  to  the  line  AD  of  the  previous 
section.   For  convenience  ARMB  is  taken  as  the  rectangle  bounded  by  s  =  s 
and  s,   and  T  =  T  and  T,   in  the  s-T  plane.   In  both  cases  the  incident  shock 
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Figure  3-^-a 
Regular  Reflection 


Figure  3.^b 
Mach  Reflection 


is  assumed  to  be  known  and  to  be  straight  so  that  the  unmodified  R-H  equations 
(l.2.l)  -  (1.2.4)  determine  the  flow  in  the  region  between  the  incident  and 
reflected  shock  regions.   RMYX  is  the  slip  stream  region.   AIDRX  and  BCMY  are 
stream  lines.   We  will  define  the  thickness  of  the  shock  regions  as 
^  T  =  length  along  the  stream  line  AI 


£\         =  length  along  the  stream  line  DR 
K 
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£s  .,     -     length  along  the  stream  line  CM  . 
M 

In  addition  let  us  define  the  lengths  along  the  arcs  RM  and  BM  as   D  and  A 

respectively.   RM  is  defined  as  the  boundary  of  the  intersection  region  such 

a  v 
that  to  the  right  of  RM  viscous  forces  such  as  w  <r=  can  be  ignored  as  they 

were  ignored  in  Section  3«^«  We  will  define  \  as  1  along  PRMQ  and  take  this 

line  as  T  =  T  . 

First  the  three-shock  problem  will  be  discussed  and  then  the  modifica- 
tions for  the  two- shock  problem  indicated,  Using  the  linear  approximation  we 
can  calculate  by  the  techniques  of  Section  3«3  the  flow  behind  the  reflected 
and  Mach  shocks  from  the  flow  ahead  of  each  and  the  shock  angles  and  curvatures 
a  ;  a  ,    k  and  k  respectively.  We  will  also  use  a  linear  approximation  of 
the  form  (3-3-1)  for  the  flow  in  the  intersection  region. 

The  method  will  be  to  determine  the  flow  along  RM  by  application  of 
the  conservation  equations  to  the  intersection  region.   Then  the  flow  is  known 
above  and  below  the  slip  stream  region  and  along  the  boundary  RM,  so  that 
the  techniques  of  Section  ^.k   will  give  four  relations  between  these  flows 
from  which  a  ,   a.  k  and  k^  are  to  be  evaluated. 

Equations  (3-3.2)  and  (3-3.3)  applied  at   s  =  s_   and  s  =  s   will 


express  a  and  A.  at  A  and  B  in  the  intersection  region  in  terms  of  5  throughout 


the  intersection  region  and  a  and  QL  j    c  _  is  known  at  the  four  corners  of  the 


intersection  region  in  terms  of  a  ,   ol^,    kR  and  k„;  therefore   its  linear 

approximation  is  known  in  ARMB  so  that  a   and  \   are  determined  in  ARMB  in  terms 

of  these  variables  and  a     and  (X..      The  energy  equation  (3.2.3)  may  be  applied 

at  s  =  s0  and  s  =  s.   to  give  p  at  R  and  M  inside  the  intersection  region 

in  terms  of  these  variables.   Then  the  momentum  equations  (3*1.2)  must  be 

applied  twice  around  suitable  contours  to  give  four  equations  that  will 

determine  A.  D,  ql,  and  a.,  in  terms  of  a„.  a...    k„  and  k...   In  the  non-viscous 

0      1  R'   M'   R      M 
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regions  the  flow  along  RX  and  MY  can  be  expressed  in  terms  of  the  flow  along  RP 

and  MQ  by  equations  (1.2.11),  and  hence  in  terms  of  a„.  a,.,  k^  and  k... 
17  '  W      M   R      M 

Finally  the  four  conditions  across  the  slip  stream  region  will  give  four 
equations  in  a  ,   0L^}    k  and  k„. 

We  will  derive  the  equations  for  a  "weak"  intersection  region,  which 
is  such  that  the  flow  variables  behave  as  indicated  for  a  weak  slip  stream  region, 
From  (3.3.2)  and  (3.3-3)  and  the  assumptions  of  the  orders  of  the  variables  we 
have 
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(3.5.1) 


±   +  0(£5) 


aA  =  a 


%    =    ao  +  0(e  } 


(3.5.2) 


We  apply  the  conservation  of  momentum  equations  (3.1.2)  separately  around  two 


halves  of  the  intersection  region  divided  by  the  stream  line  s  = 


so  !  si 


To  take  account  of  the  flow  above  s  =  s   and  below  s  =  sn  we  will  evaluate 
the  integrals  along  these  lines  in  the  regions  adjacent  to  the  intersection 
region,  namely  the  shock  and  non- viscous  regions.   The  viscous  terms  can  once 
again  be  neglected.  We  get  by  using  the  trapezoidal  rule  for  integration 
around  the  upper  half  of  ARMB 
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3E 


where  n  =  3  and  n  =  2,       C  x\  means  X.  +  3X  ,  and  Y.  =  pu.  — -;  and 


around  the  lower  half  of  ARMB 
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Taking  three  times  the  (3° 5° 3)  minus  (3.5.^)  we  get 
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and  taking  three  times  (3.5A)  minus  (3« 5 <•  3 )  we  get 
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(3.5.6) 


We  now  evaluate  the  terms  on  the  left  hand  sides  to  get  the  four 
equations 
2Dt*(  \2R  -  \      )   ♦  (pR  -  p  )rt  ]  +  SAtSPj  "  PA  -  PR  -  PMJ 


■  "3\(PR  "  PI)  +  3,iI(pI  "  PA>  -AM(pM  "  PA>  +  °{i     >  (3-5-7) 
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"  VPR  "   PI>  "  V*I  -   PA>   +  3VPM  '   Pa)   +  °(*     '  (3-5-S) 


*+  1     _,  w   PR  '   PA 

— -  D(pR  -   pA)(  — — ■ 


lh   >+i*(W+  WW    3Ul> 


-    3V|2r"  W"  3V^t      >pa)+V^--   ^.)  +  °(^3) 


'21         '2A'  MV?2M        7  2A' 


(3.5-9) 


and 


Pw     ~      Pa  ^2 

7-     D(%  -   PA)(  — y yTl    >  +  2A(§2I  +    &2A  •     52R  "    W 


i.( 


»+AI^2I-     ^2A»-    3Vl 


Rv  >2R         721'  Ix  721         72A 


2M         ?2A 


)   +   0(£3) 


(3.5.10) 
We  have  used  the  assumption  that  the  pressures  p  evaluated  at  the  same  point  in 
either  the  slip  stream  region  or  the  non- viscous  region  agree  to  terms  in  &     . 
This  can  be  verified  from  the  energy  equation  (3.2.3)  which  implies  that 
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for  the  stream  line  s  =  sn  or  s   considered  as  belonging  to  either  the  inter- 
section region  or  the  shock  and  non- viscous  regions.   The  conservation  of  mass 


0  = 


__              _ 

T 
1 

J     n_ 

= 

r 

+  o(*4) 


(3.5.12) 


'0 


is  also  satisfied  on  s  =  sn  or  s   in  both  regions.   Since  p  ,    v  and  p  at 


T  =  T   are  independent  of  the  region  of  evaluation  for  s  =  s_  or  s,  and  the 
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80 
same  is  true  for  v  at  T  =  T  ,   (3.5.11)  and  (3. 5. 12)  imply  that  p  and  P  are 
independent  of  the  region  of  evaluation  at  R  or  M  except  for  terms  of  0(fe  J 
and  higher. 

3.6       In  this  section  we  will  apply  the  modified  and  unmodified  R-H  equations 
and  the  modified  matching  conditions  to  determine  the  shock  region  pattern. 
Let  the  oncoming  flow  be  scaled  so  that 

*o  ■  Po  "  1 
1  A2 

assume  that  the  angle  of  the  incident  shock  is 

if 

a  =    —   +  afc 

and  its  curvature  zero.  The  flow  behind  the  incident  shock  region  is  given  by 
the  R-H  equations  (l.2.l)  -  (1.2.U).  The  Mach  number  behind  the  shock  is 

m2  =  1  +  (2a2  -  l)£2 
J. 

and  the  direction  of  the  flow  is 


^      "*•"         2     f.  2x,3 


21    2        Tf+  1  a(l  "  a  ^ 
radians;  to  the  horizontal.   Let  us  define  the  reflected  shock  angle  by 

2      2       2 
where  6R  =  (2a  -  l)6    and  the  Mach  shock  angle  by  a^  =<T/2  +   0  £  .  Then  the 

modified  R-H  equations  (3-3.17)  -  (3.3.2*0  define  the  flow  behind  the  shock 

regions  in  terms  of  S.  k,  6  and  k^  with  suitable  substitutions.   In  the  non- 

M         R 

viscous  regions  we  can  write 
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Tt         ij   a>  s 


from  (1.2.11).   Let  a  =0,   a_=  3,  6M  =6;   then  for  Z  =  R  or  M  we  have 
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Equations  (3.5-9)  and  (3-5-10)  are  linear  equations  in  D  and  A  ,  and  thus 
express  D  and  A  in  terms  of  (3,  o  ,  k  and  k^.   Substituting  these  relations 
into  (3 • 5 • 7 )  and  (3-5.8)  we  get  b  and  b  in  terms  of  these  variables.   Sub- 
stituting these,  (3.6.I)  and  (3-6„2)  into  (3-^-iM  and  (3-^-5)  we  get  two  rela- 
tions, and  substituting  (3-6.1)  and  (3-6.2)  into  (3.^-6)  and  (3-^-7)  we  get  two 
more  relations  in  Pj  o  ,  k  and  k^.  Thus  we  have  four  equations  to  solve  for 
these  four  variables.   One  solution  that  can  be  seen  immediately  is  k  =  1c.  =  0, 
3  and  6  given  by  the  solution  of  the  non- viscous  equations.   This  is  obvious 
since  k  =  k^  =  0  implies  that  derivatives  with  respect  to  T  are  0  behind  the 
shocks,  and  hence  (3-^-6)  and  (3-^-7)  are  satisfied,  leaving  (3-^«^-)  and  (3-^-5) 
for  which  the  right  hand  sides  vanish  so  that  they  give 
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=  0 


which  are  the  usual  non- viscous  conditions.   Solutions,  if  they  exist,  for  the 

variables  (3,  o  ,  kR,  k^,  D,  &  ,  b  and  b  have  the  order  that  was  assumed,  so 

that  the  assumptions  are  consistent.  The  thicknesses  of  the  shocks  A „■ A.,  and 

A  are  proportional  to  w,  so  from  equations  (3-5-7)  -  (3-5-10)  we  have  that  A 


82 
and  D  are  proportional  to  w  while  b  and  b  are  independent  of  it.   Then  from 
(3.^.1j.)  _  (3.^.7)  we  find  that  k  and  k^  are  inversely  proportional  to  w  hut 
P  and  o  are.  independent  of  it.   If  a  solution  of  these  equations  exists  with 
either  k  or  k^  non-zero,  then  as  the  viscosity  w  tends  to  zero,  the  non-zero 
k's  will  tend  to  co  „  Therefore  we  cannot  have  a  singularity  of  this  type  in 
which  the  non-viscous  R-H  equations  and  matching  conditions  break  down  unless 
k  is  infinite  when  the  viscosity  is  zero. 

3.7       For  regular  reflection  in  viscous  flow,  there  are  two  cases;  reflection 
at  a  bounding  wall  and  reflection  at  an  axis  of  symmetry.   In  non-viscous  flow 
these  cases  are  identical.   In  reflection  at  a  bounding  wall  in  viscous  flow 
the  velocity  is  prescribed  as  zero  on  the  wall.   Since  the  fluid  away  from  the 

wall  is  moving  there  will  be  a  boundary  layer  acrosss  which  the  velocity  is  non- 

[5] 
linear  (see  Goldstein   ).   The  reason  for  the  non-linearity  is  that  at  the  wall 

b  v 

the  w  \- —  must  be  large  so  that  viscous  drag  can  slow  the  fluid  while  at  the 

o  v 
non- viscous  boundary  of  the  boundary  layer,  w  -v—  must  be  small.  Therefore, 

our  linear  approximation  to  v  is  inappropriate  to  reflection  at  a  bounding 

wall.   We  will  only  consider  reflection  at  a  line  of  symmetry  for  which  the 

only  condition  is  that  r   =  7T/2  along  the  stream  line  MY. 

We  replace  the  variables  3  and  k  by  the  velocity  and  its  T  derivative 

f^2\    "£    tr/ 
at  M  and  write  I  -5-=    =  ?       -  '1/2   =  0.   Using  these  in  the  three  shock  equa- 
tions and  writing  &  =A  we  will  get  a  similar  set  of  equations  appropriate  to 
regular  reflection.  All  remarks  concerning  the  behavior  of  Q,  k_,  &  ,  D, 

b  and  b  in  the  three  shock  problem  are  applicable  to  these  variables  in  the 

[9] 
two  shock  problem.   Since  Smith    has  observed  regular  reflection  for  weak 

shocks  at  incident  angles  for  which  no  solution  of  the  R-H  equations  exists, 

we  hope  that  infinite  curvature  of  the  shock  will  allow  solutions  to  these 

equations  to  exist  when  the  viscosity  tends  to  zero. 
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IV  CONCLUSION 

In  Chapter  II  we  studied  the  problem  of  a  wedge  in  an  oncoming  super- 
sonic non- viscous  flow  such  that  the  shock  was  attached  to  the  wedge.   This  was 
similar  to  the  problem  of  regular  reflection.  We  saw  how  by  the  method  of 
characteristics  we  could  solve  for  the  flow  behind  the  shock  when  the  flow  ahead 
was  uniform,  the  wedge  wall  straight  and  the  flow  initially  supersonic.   Similar 
methods  will  suffice  to  treat  the  problem  when  the  flow  ahead  is  non-uniform  and 
the  wedge  wall  curved.   When  the  flow  emerging  from  the  shock  was  initially  sub- 
sonic an  analytic  solution  existed  only  for  a  countable  number  of  wedge  angles. 
For  other  angles  a  class  of  singular  solutions  of  the  differential  equations  was 
found  that  satisfied  the  local  boundary  conditions,  and  this  class  was  seen  to 
have  its  analogy  in  pseudo- stationary  flow.   One  problem  that  was  not  discussed 
was  the  use  of  other  boundary  conditions  to  determine  the  parameters  of  the 
solution.   Until  this  has  been  treated,  no  discussion  of  uniqueness  is  possible. 

The  solution  had  a  branch  point  singularity  at  the  origin,  similar  to 
the  singularity  used  by  Guderly    (Chapter  VI )  in  his  discussion  of  an 
irrotational  flow  behind  a  shock.   A  branch  point  is  also  found  in  two-dimensional 
potential  problems  at  an  angle  in  an  e qui- potential  line.   The  similarity  is  due 
to  the  fact  that  the  nonlinear  equations  of  general  or  irrotational  flow  can  be 
linearized  for  small  deviations  of  the  dependent  variables  around  the  origin. 
Since  the  equations  are  elliptic,  a  suitable  linear  transformation  of  the  inde- 
pendent variables  will  then  transform  them  into  the  Laplace  equation. 

For  some  incident  shocks,  regular  reflection  was  no  longer  possible. 
This  corresponded  to  the  wedge  angle  being  too  large  for  a  shock  to  be  attached 
to  it.   Experimentally,  either  regular  reflection  or  a  three  shock  configuration 
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with  a  slip  stream,  called  Mach  reflection,  is  observed.  For  weak  incident 
shocks,  no  solution  of  the  algebraic  R-H  equations  and  matching  conditions 
across  the  slip  stream  existed  so  we  could  not  use  the  techniques  of  Chapter  II 
which  assumed  the  existence  of  such  solutions.   Consequently  we  examined  the 
equations  of  viscous  flow.  The  problem  was  simplified  by  replacing  the  mathe- 
matical discontinuities,  called  shocks  and  slip  streams,  of  non-viscous  flow 
by  regions  in  which  viscosity  played  a  role.   In  other  regions  the  viscosity 
was  neglected. 

We  first  examined  the  structure  of  a  plane  two-dimensional  shock 
region  using  a  linear  approximation.  We  derived  the  equations  when  the  Mach 
number  of  the  flow  ahead  was  close  to  1.  We  are  particularly  interested  in 
this  case  because  it  corresponds  to  those  problems  where  the  largest  divergence 
between  theory  and  experiment  has  been  observed.   Next  we  examined  a  viscous 
region  corresponding  to  a  slip  stream  in  a  similar  manner, and  then  applied  the 
equations  derived  to  the  problems  of  two  and  three  shock  intersections.   The 
equations  obtained  admitted  the  usual  non-viscous  solutions  when  the  shock 
curvature  was  zero.   Newton's  method  in  the  four  variables  (3/0,  k„j  and  k  , 
was  used  to  study  the  equations  numerically.   The  only  solution  found  prescribed 
a  negative  shock  in  violation  of  the  second  law  of  thermodynamics. 

If  these  equations  have  no  physically  acceptable  solutions  we  must 
re-examine  the  approximations.  The  first  step  should  be  to  use  quadratic 
approximations  in  s  for  the  slip  stream  and  intersection  regions  since  we  saw 
that  a  linear  approximation  was  inappropriate  to  the  study  of  boundary  layers . 

In  the  solutions  of  the  linear  approximations,  if  they  exist,  the 
variables  k,  the  curvatures  of  the  shocks,  enter  only  in  the  form  kw.   In 
Chapter  II  we  discussed  the  solution  of  pseudo- stationary  flow  equations  where 
the  variables  of  length  x  and  x  are  replaced  by  the  variables  a  =  x  /t  and 
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a  =  x  /t  where  t  is  the  time.   The  pseudo- stationary  shock  curvature  k  =  -r —  , 
where  s  has  the  dimensions  of  length/time.   Therefore,  terms  in  kw  will  be 

functions  of  a  ,    a  and  t. 

[9] 

This  is  in  contradiction  to  the  experiments  of  Smith    and  others , 

who  have  observed  no  measurable  deviation  from  pseudo- stationary  behavior  in 
their  shock  tube  experiments.   Thus  suggests  that  the  region  over  which  viscosity 
plays  a  role  is  too  small  to  be  observed  in  their  experiments  or  that  there  does 
exist  a  form  of  the  solution  in  which  the  region  affected  by  viscosity  depends  on 
a  and  a  only. 

A  method  of  applying  the  approximations  for  the  shock  regions  con- 
sistently to  a  boundary  value  problem  is  needed  before  we  can  come  to  understand 
the  effect  of  the  bounding  walls  on  the  flow  patterns  for  various  viscosities. 
In  view  of  the  fact  that  the  equations  are  different  from  the  non- viscous  theory, 
the  surprising  result  is  that  there  does  exist  a  solution  corresponding  to  the 
non-viscous  theory. 
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APPENDIX  I 


The  viscous  term  in  the  energy  equation  across  a  weak  shock  region  is 
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where  the  terms  are  evaluated  somewhere  in  (Tn,  T  )  unless  otherwise  indicated, 
However  "by  our  assumptions 
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therefore  the  viscous  term  has  order 
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A  similar  proof  shows  that 
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